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PREFACE 

This  report  summarizes  the  activities  supported  under  the  Office  of  Naval  Research  (ONR) 
Contract  No.  N00014-89-C-0117.  The  ONR  technical  monitor  for  the  activity  was  Dr.  Frank 
Herr.  The  ERIM  principal  investigator  was  Dr.  Robert  Shuchman.  Dr.  Robert  Onstott  and  Mr. 
Scott  Gaboury  collected  the  active  and  passive  microwave  data  during  SLICKEX,  while  Dr. 
Clifford  Rufenach  and  Dr.  Robert  Shuchman  aided  Dr.  Robert  Onstott  in  the  analysis  and 
interpretation  of  the  data. 
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1.0  EXECUTIVE  SUMMARY 

We  measured  the  radar  backscattered  cross  section  across  naturally  occurring 
ocean  slicks  during  SLIX  ’89  and  also  developed  expressions  for  the  sensitivity  of 
active/passive  sensors  for  extracting  sea  surface  temperature. 

The  results  include  ship-based  radar  scatterometer  measurements  at  L-band  and  X- 
band  and  microlayer  chlorophyll  fluorescence  measurements  on  the  southern 
California  continental  shelf  across  several  surface  slicks  during  October,  1989.  These 
measurements  allowed  the  simultaneous  extraction  of  both  the  short  wave  damping  at 
Bragg  wavelengths  of  24  cm  and  4  cm  and  the  chlorophyll  A  biological  concentration 
with  sufficient  spatial  resolution  to  observe  features  within  individual  slick  bands. 

The  measurements  show:  (1)  the  largest  peaks  in  fluorescence,  a  measure  of  the 
peak  Chlorophyll  A  concentration,  are  co-located  with  the  short  wave  damping  at  24 
cm  and  4  cm  wavelengths,  (2)  the  maximum  damping  ratio,  between  the  24  cm  and  4 
cm  waves  is  at  least  16  for  wind  speeds,  3-4  m/s,  and  (3)  the  radar  roughness 
signatures  at  these  short  Bragg  gravity  wavelengths  suggest  the  slicks  are  caused  by 
damping  rather  than  straining  of  these  waves.  Additional  shipboard  measurements  are 
needed  to  further  quantify  these  preliminary  results. 

We  study  the  feasibility  of  using  a  dual  linear-polarized  radiometer  operating  at  6 
GHz  in  conjunction  with  a  scatterometer  as  a  method  c .  extracting  sea  surface 
temperature.  A  model  originally  developed  by  Wilheit  [1979]  is  used  to  investigate 
the  sensitivity  of  radiometric  temperature  to  changes  in  wind  speed,  U,  at  a  height  of 
20  m;  that  is,  the  sea  temperature  accuracy  depends  directly  on  the  radiometric 
temperature  accuracy.  These  changes  are  corrected  using  this  extracted  wind  speed 
based  on  a  scatterometer  wind  accuracy  of  ±  1  m/s.  We  show,  as  previously 
reported,  that  a  vertically  polarized  radiometer  is  less  sensitive  to  wind  speed  changes 
than  a  horizontally  polarized  one  near  an  incidence  angle  0  =  60°.  Furthermore,  we 
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quantify  the  wind  variation;  ATb^AU  <  0. 12°C/m/s  when  U  >  10  m/s  and 
ATb”/AU  <  2.0°C/m/s  when  U  >  6  m/s  for  vertical  and  horizontal  polarization, 
respectively.  Thus,  the  resulting  error  is  ±  0. 12°C  and  ±  2.0°C  for  Tb'"  and  Tg” 
based  on  the  extracted  wind  accuracy  ±  1  m/s.  These  results  suggest  that  a 
scatterometer  can  be  used  to  correct  the  wind  speed  dependence  of  a  vertically 
polarized  radiometer  operating  near  6  GHz  using  a  pencil  antenna  beam  pointing  at  an 
incidence  angle  near  55°. 

A  number  of  symposium  and  referred  journal  articles  have  resulted  from  this 
research  effort.  Table  1  lists  these  articles,  while  the  appendix  of  this  report  provide 
a  copy  of  each. 
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2.0.  INTRODUCTION 

Numerous  active  (scatterometric)  and  passive  (radiometric)  microwave 
measurements  of  the  ocean  surface  have  been  acquired,  analyzed  and  interpreted  over  the 
last  twenty  years.  However,  few  simultaneous  active  and  passive  measurements  have 
been  interpreted  in  terms  of  physical  properties  of  the  ocean  surface.  Simultaneous 
measurements  have  the  potential  of  extracting  sea  surface  temperature  more  accurately 
than  passive  measurements  alone.  Furthermore,  simultaneous  shipboard  microwave 
measurements  in  conjunction  with  the  Self  Contained  Underway  Microlayer  Sampler 
(SCUMS)  towed  from  the  side  of  a  ship  [Carlson  et  al.,  1988]  during  the  Marine 
Microlayer  Slick  Experiment  (SLIX  ’89)  produced  information  about  the  concentration 
of  the  marine  microlayer  material  and  associated  damping  of  the  short  gravity  waves 
within  individual  ocean  slick  bands. 

The  radiometric  temperature  of  an  ocean  surface  depends  primarily  upon  its  sea 
surface  temperature,  salinity,  surface  wind  speed  (surface  roughness),  and  foam  cover. 
A  smooth  surface  does  not  depend  on  the  wind  speed  or  foam  cover;  indeed,  foam  cover 
is  not  important  for  wind  speeds  less  than  about  7  m/s.  A  smooth  sea  surface  depends 
on  surface  temperature  and  salinity.  This  dependence  is  known  accurately  through  its 
relationship  with  the  complex  dielectric  constant  of  sea  water  [Klein  and  Swift,  1977]. 
The  surface  wind  influences  the  radiometric  temperature  through  the  ocean  wave 
spectrum.  However,  measurements  are  not  fully  consistent  with  each  other  or  with  the 
theory  especially  when  foam  is  present  [Hollinger,  1971  and  Swift,  1974].  Additional 
surface  based  measurements  are  required  to  validate  radiometric  models  before  additional 
aircraft  and/or  spacecraft  based  measurements  can  be  fully  interpreted. 

Radar  scatterometers  measure  the  ocean  surface  roughness  at  Bragg  wavelengths, 
3  cm  -  30  cm,  within  the  radar’s  antenna  beamwidth.  The  SCUMS  platform  contains 
Fluorescence  and  UV  absorbance  detectors.  Fluorescence  detectors  measure  the 
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concentration  of  the  biological  material,  chlorophyll  A;  whereas,  the  UV  absorbance 
detectors  (280  nm)  measure  the  concentration  of  chemical  material.  The  ability  of  the 
microlayer  to  dampen  short  waves  increases  with  the  concentration  of  the  microlayer 
material.  This  concentration  varies  across  slicks  due  to  the  horizontal  convergence  of 
ocean  currents  at  the  surface  and/or  by  horizontal  convergence  due  to  wind  stress.  That 
is,  the  surface  roughness  near  and  within  slick  bands  vary  with  ocean  current  and  wind. 

Slicks  on  the  southern  California  and  Baja  continental  shelf  can  be  due  to  internal 
waves  currents,  small  scale  ocean  circulation  and  wind  stress.  The  internal  waves 
typically  travel  along  the  thermocline  at  the  same  speed  as  the  companion  slick  bands. 
These  slicks  are  characterized  by  a  series  of  slick  packets;  each  packet  separated  by  a 
distance  of  15  -  20  km  with  an  envelope  of  quasi-period  parallel  slick  bands  with 
wavelengths  of  about  3(X)  m  [Ewing,  1950;  and  Apel  and  Gonzalez,  1983]. 

The  slick  characteristics  can  be  used  to  distinguish  internal  wave  generated  slicks 
from  other  slick  generating  mechanisms.  The  surface  roughness  at  short  gravity 
wavelengths  associated  with  quasi-periodic  slicks  have  been  observed  by  a  number  of 
workers  using  different  measuring  methods;  see,  e.  g.,  Ewing  [1950],  LaFond  [1962], 
Rufenach  and  Smith  [1985],  Alpers  [1985],  Shuchman  et  al.  [1988],  Lyzenga  and  Benr.ett 
[1988]  and  Gasparovic  et  al,  [1988].  Alternating  bright  and  dark  bands  in  radar  images 
have  been  interpreted  as  horizontal  current  gradients  produced  by  internal  waves;  see, 
e.  g.,  Alpers  [1985];  or  equivalently,  increases  and  decreases  in  radar  cross  section 
(RCS).  In  contrast,  RCS  decreases  are  directly  related  to  the  damping  of  these  short 
waves  due  to  increases  in  the  concentration  of  the  surface  marine  microlayer  material. 
Surface  roughness  changes  have  usually  been  interpreted  in  terms  of  either  damping  or 
straining  of  the  short  waves  but  not  both.  We  believe  that  both  mechanisms  may  be 
important  on  some  occasions. 

The  concentration  of  the  surface  microlayer  material  has  been  related  to  the  bulk 
viscosity  [Carlson,  1982;  Carlson,  1987].  Microlayer  viscous  damping  is  caused  by  either 
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the  bulk  viscosity,  syrupy  nature  of  the  surface  layer,  or  the  elasticity,  visco-elastic 
nature  of  this  thin  surface  layer.  The  elastic  properties  are  equivalent  to  an  additional 
viscosity  of  the  surface  microlayer  [Levich,  1962].  We  do  not  attempt  to  distinguish 
between  these  two  viscous  damping  mechanisms. 

Detection  of  bands  of  slicks  associated  with  internal  waves  by  passive  microwave 
sensors  have  not  been  reported.  Passive  microwav'  sensors  have  been  primarily  operated 
fiom  spacebome  platforms  with  a  spatial  resolution  >  20  km.  This  resolution  could 
explain  the  lack  of  detection  since  the  slick  packet  spacing  is  typically  k  -  15  km. 
Increases  in  emissivity  have  been  measured  across  man  made  oil  slicks  by  microwave 
sensors  [Hollinger  and  Mannella,  1973]  whereas  one  example  of  a  large  emissitivity 
decrease  has  been  reported  across  a  manmade  monomolecular  slick  [Huhnerfuss  and 
Alpers,  1983].  Additional  shipboard  multi-frequency  active  and  passive  microwave 
measurements  in  conjunction  with  in  situ  measurements  across  manmade  and  natural 
slicks  are  needed  to  better  understand  the  radio  perturbations  caused  by  surface  slicks 

This  report  is  limited  to  determining  an  optimum  radiometric  frequency  for  the 
measurement  of  sea  surface  temperature  using  an  active/passive  sensor  system  and  the 
comparison  of  scatterometer  and  radiometer  measurements  with  in  situ  fluorescence 
measurements  near  two  surface  slick  groups  off  the  California  coast. 
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3.0.  SLIX  ’89  EXPERIMENTAL  DESCRIPTION 

The  ship-board  SLIX  ’89  designed  to  traverse  surface  slicks  off  the  coast  of 
southern  California  near  San  Diego  during  a  two  week  period  in  October  1989.  The  ship 
heading  during  the  experiment  was  selected  to  intersect  slick  groups  visually  observed 
from  the  bridge  of  the  ship.  This  experiment  included  a  variety  of  in  situ  and  remote 
sensing  instruments  including  the  microwave  sensors  and  SCUMS  catamaran  in  situ 
detectors. 

The  microwave  sensors  were  mounted  about  10  m  above  the  sea  surface  toward 
the  bow  of  the  ship,  R/V  Wecoma,  whereas  the  in  situ  instruments  were  towed  by  a 
catamaran  attached  to  a  boom  on  the  port  side  of  the  ship.  The  wind  velocity  at  a  height 
of  about  10  m  was  measured  near  the  bow  of  the  ship  and  sea  surface  temperature  was 
measured  at  a  depth  of  15  cm.  We  operated  the  following  microwave  sensors: 

(1)  a  dual  wavelength  radar  scatterometer,  and 

(2)  a  dual  wavelength  radiometer, 

both  mounted  on  the  port  side  with  their  antennas  pointing  perpendicular  to  the  ship’s 
track.  The  radar  operating  wavelengths  were  L-band  and  X-band  (20  cm  and  3.2  cm) 
while  the  radiometers  operated  at  C-band  and  Ka-band  wavelengths  (5.5  cm  and  0.86 
cm).  The  radiometer  outputs  were  carefully  monitored  for  interference  from  the 
scatteiometer.  The  relevant  instrument  parameters  are  given  in  Table  1.  The  radiometers 
were  calibrated  using  a  hot  absorber  and  cold  (near  zenith  sky  scan)  load. 

The  in  situ  measurements  of  interest  were  taken  on  board  the  catamaran  towed 
about  20  m  aft  of  the  microwave  sensors  from  the  port  side  of  the  ship,  SCUMS  [Carlson 
et  al.,  1988].  The  UV  absorbance  is  related  to  the  concentration  of  the  absorbing  organic 


6 


?>ERIM 


TABLE  1. 

Scatterometer  and  Radiometer  System  Parameters 


RADAR  RADIOMETER 

FREQUENCY 

1.5  GHz 

9.38  GHz 

5.4  GHz 

35  GHZ 

WAVELENGTH 

20  cm 

3.2  cm 

5.5  cm 

0.86  cm 

ANTENNA  SPOT  SIZE 
at  0  =  55° 
at  0  =  25°  § 

3.7  m 

1.5  m 

1.3  m 

0.5  m 

m 

m 

INTEGRATION  TIME 

50  sec 

50  sec 

50  sec 

50  sec 

BANDWIDTH 

350  MHz 

575  MHz 

70  MHz 

500  MHz 

POLARIZATION 

HH,VV  or  VH 

HH.VV  or  VH 

H  or  V 

H  or  V 

RESOLUTION 

1.1  dB 

2.6  dB 

0.21  “K 

0.12°K 

ABSOLUTE 

ACCURACY 

±  1  dB 

±  1  dB 

±  1°K 

±  TK 

*  H  signifies  horizontal  polarization  and  V  signifies  vertical  polarization. 

§  incidence  angle  used  during  measurements  across  slick  groups  "A"  and  "B" 
on  October  16. 
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molecules,  whereas  the  chlorophyll  A  fluorescence  is  related  to  the  concentration  of  the 
biota.  The  SCUMS  measurements  include:  the  fluorescence  and  the  UV  absorption 
measured  in  both  the  microlayer  and  sub-layer.  The  microlayer  material  was  attached  to 
a  rotating  drum  then  skimmed  from  the  drum  surface  whereas  the  sub-surface  material 
was  collected  at  a  depth  of  about  15  cm.  The  microlayer  is  thought  to  be  50-100  ^m 
thick.  The  fluorescence  and  UV  absorbance  show  similar  variations  across  the  slicks 
implying  that  the  concentration  of  the  surface  microlayer  material  can  be  represented  to 
first  order  by  either  detector. 

A  side  view  of  the  ship,  catamaran  and  remote  sensing  instruments  are  shown  in 
Figure  1(a)  and  a  side  view  of  the  catamaran  and  associated  instrumentation  is  shown  in 
Figure  1(b).  A  photograph  of  the  microwave  antennas  are  given  in  Figure  2.  Selected 
measurements  during  SLIX  ’89  are  given  in  the  following  section. 

3.1  MEASUREMENTS 

Two  surface  slick  groups  off  the  coast  of  San  Diego  during  October  16  were 
selected  for  analysis.  These  two  groups  each  exhibit  three  slick  bands  during  a  period 
when  both  the  microwave  and  in  situ  instruments  were  operating.  Furthermore,  these  two 
slick  groups  were  visually  observed  from  the  ship  deck  as  recorded  in  the  operational 
log.  We  distinguish  between  groups  and  packets  in  the  present  work:  packets  are 
associated  with  internal  waves;  typically  separated  by  15  -  20  km.  Groups  are  more 
closely  spaced  clusters  of  slick  bands,  several  occurring  within  a  packet,  due  to  the 
complex  nature  of  the  bands  or  generating  mechanisms  other  than  internal  waves.  The 
orientation  of  these  slick  bands  relative  to  the  radar  look  direction  is  not  known  since 
radar  imagery  was  not  available  and  no  orientation  information  was  recorded  in  the  log. 
However,  the  experiment  objective  was  to  traverse  the  slick  bands.  That  is,  the 
fluorescence  and  radar  measurements  show  a  fluctuating  component  along  the  ship  track. 
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91-20063 

Figure  1.  (a)  Artists  Rendition  of  the  Research  Vessel  Wecoma  Relative  to  Towed 

Catamaran  (SCUMS).  The  Microwave  Equipment  Measures  Radar  Cross 
Section  (RCS)  and  Radiometric  Brightness  Temperature.  The  SCUMS 
Instrument  Measures  Fluorescence  and  UV  Absorbance  of  the  Surface 
Layer  Organic  Material,  (b)  Side  View  of  the  Catamaran  and  Attached 
SCUMS  Sampler.  The  Micro-Layer  Film  Attaches  Itself  to  the  Rotating 
Cylindrical  Mirror  Surfaced  Drum.  The  Microlayer  Material  is  Extracted 
From  the  Drum  and  Pumped  Into  the  Ship  Where  it  is  Analyzed. 
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Figure  2.  Photograph  of  the  Radio  Microwave  Sensors:  Radar  Scatterometers 
Passive  Radiometers. 
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We  have  no  direct  measurements  in  the  water  column  that  these  slicks  were 
produced  by  internal  waves  or  any  other  mechanism.  However,  the  separation  between 
the  two  groups  (*  30  min.)  and  a  slick  band  spacing  (3  -  5  min.),  see  Figure  3,  suggest 
spatial  separations;  the  order  of  a  few  km  and  hundreds  of  meters,  respectively.  The  3  - 
5  min.  quasi-periodic  spacing  are  consistent  with  internal  wave  produced  slick  bands 
provided  the  ship  speed  (0.5  -  1  m/s)  is  much  larger  than  that  component  of  the  slick 
velocity  along  the  ship  track.  The  spatial  separation  of  a  few  kilometers  suggest  that  the 
two  slick  groups  are  not  internal  wave  associated  packets.  The  grey  areas  depicted  in 
Figure  4  give  the  location  of  the  two  slick  groups  relative  to  the  ship  track  during 
October  16.  The  two  slick  groups  occurred  in  the  morning  hours  between  10:40  to  11:20 
Pacific  Daylight  Time  (PDT).  There  were  no  clouds  or  rain  during  this  period. 

An  overview  of  the  microwave  and  in  situ  measurements  using  stacked  plots  of 
the  various  parameters  is  given  in  Figure  4.  The  measured  parameters  (1  min  averages) 
are  displayed  on  the  ordinate  and  the  time  in  minutes  is  displayed  on  the  abscissa  in  all 
the  plots.  The  values  of  RCS  (o"),  fluorescence  relative  to  the  sub-layer  (F,)  and 
radiometric  brightness  temperature  (Tb)  were  selected  for  display.  The  two  o"’s  at  L-band 
and  X-band  are  indicated  in  Figure  4(a).  The  relative  slick  fluorescence  (F,)  is  obtained 
by  dividing  the  microlayer  fluorescence  by  the  sub-layer  fluorescence.  The  inherent  delay 
in  skimming  and  pumping  the  surfactant  microlayer  has  been  corrected.  The  delay  due 
to  the  spatial  separation  of  the  microwave  instruments  and  the  catamaran,  ^  20  sec,  has 
not  been  corrected.  The  radiometric  brightness  temperatures  at  C-band  and  Ka-band  show 
no  discernible  change  across  the  slick  bands.  That  is,  any  surface  slope  decrease  due  to 
the  slicks  was  not  detected  by  these  radiometers. 

The  UV  absorbance  and  the  fluorescence  show  similar  variations  across  the  two 
surface  slick  groups  of  October  16.  We  select  the  fluorescence  for  detailed  analysis  in 
the  present  work.  The  fluorescence  changes  from  10:40  to  11:20  is  a  measure  of  the 
concentration  of  the  chlorophyll  A  within  the  two  slick  groups  separated  by  a  relatively 
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Pacific  Daylight  Time  (PDT) 

91-2CX)67R1 

Figure  3.  Stacked  Plots  of  Measured  Parameters  for  the  40  Minute  Period  10:40  to 
11:20  PDT  for  Two  Slick  Wave  Groups  From  R/V  Wecoma,  Morning  of 
October  16,  1989.  The  Measured  Parameters  are:  (a)  Radar  Backscatter 
Cross  Section  ((f),  (b)  Slick  Fluorescence  Relative  to  the  Sub-Layer  (F,)  and 
(c)  Radiometric  Brightness  Temperature  (Tg).  The  Polarization  was  Vertical 
for  the  Radiometer  and  Horizontal  Transmit  and  Horizontal  Receive  for  the 
Radar  Both  at  an  Incidence  Angle  of  0  =  25®. 
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Figure  4,  Plan  View  of  Ship  Track  and  Areas  of  Surface  Slick  Group  on  the  R/V 
Wecoma  on  October  16,  1989.  The  Ocean  Depth  is  Indicated  by  the 
Contours  90  m  and  365  m. 
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surfactant  free  region.  The  first  group  is  designated  as  group  "A"  and  the  second  group 
is  designated  as  group  "B".  The  three  peaks  in  fluorescence  F,  within  each  slick  group 
are  interpreted  as  three  slick  bands. 

The  relative  fluorescence  F,  increased  to  a  peak  value  of  7  within  the  slick  group 
"A"  near  10:47.  The  radar  cross  section  also  decreased  about  10  dB  at  X-band  and  6  dB 
at  L-band,  both  near  10:47,  co-located  with  the  fluorescence  peak.  This  co-location  is 
direct  evidence  that  the  short  Bragg  waves  are  damped  at  the  peak  surfactant 
concentration  near  the  center  of  the  slick  band. 

The  relative  fluorescence  F,  increased  to  a  peak  value  of  8  within  slick  group  "B" 
near  11:10.  Different  RCS  decreases  are  apparent  within  slick  group  "B"  compared  to 
slick  group  "A".  The  RCS  decreased  about  10  dB  at  X-band  similar  to  slick  group  "A" 
whereas  the  cross  section  decreased  about  2  dB  at  L-band  in  contrast  to  the  6  dB 
decrease  in  slick  group  "A".  Furthermore,  the  X-band  decrease  occurred  over  a  longer 
interval,  about  8  min,  than  group  "A".  This  large  damping  at  X-band,  reducing  the 
backscattered  energy  to  near  the  noise  floor  of  the  scatterometer,  does  not  allow  co- 
location  of  the  peak  fluorescence  F,  with  the  radar  minima  signatures.  That  is,  this  X- 
band  damping  occurred  across  several  slick  bands  within  Group  "B",  whereas,  the  L- 
band  radar  minima  are  approximately  co-located  with  all  three  fluorescence  peaks. 

The  wind  speed  changed  from  about  1.5  m/s  to  about  4  m/s  while  the  wind 
direction  changed  from  about  60°  to  about  30°  with  respect  to  the  radar  look  direction 
from  group  "A"  to  group  "B".  The  wind  was  blowing  generally  toward  the  south  and 
south-east  toward  the  radar  look  direction.  The  ship  speed  was  about  0.5  m/s.  The  wind 
speeds  were  light,  1.5-3  m/s  inside  group  "A"  at  a  wind  direction  of  about  60°.  Slightly 
higher  wind  speeds,  3-4  m/s  were  observed  inside  group  "B"  at  a  direction  of  about  30°. 
The  differences  between  the  X-band  RCS  decreases  within  the  two  slick  groups  could  be 
due  to  the  wind  direction  changes  between  the  two  slick  groups.  These  differences  are 
discussed  in  more  detail  in  the  following  sections.  Decreases  in  RCS  co-located  with 
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fluorescence  peaks  were  also  observed  within  several  other  slicks. 

3.2  BRAGG  SCATTERING  MODEL 

Typical  ocean  remote  sensing  radars  operate  at  microwave  wavelengths  in  the 
centimeter  or  decimeter  range  at  incidence  angles,  0  from  20-70°  within  which  the  radar 
backscatter  is  dominated  by  Bragg  scattering.  The  RCS  to  first  order  for  horizontal 
transmit  and  receive  polarization  is  given  by  [Valenzuela,  1978], 

AK)  =  4  ko^  cos^  0  I  g„H(0)  1  '  S(k„0)  (Eq.  1) 

where  ko  is  the  radar  wavenumber,  S(k„ky)  is  the  spectral  density  of  the  ocean  surface 
waves  in  rectangular  coordinates,  and  gHwC©)  is  the  first  order  scattering  coefficient 
dependent  on  the  relative  dielectric  constant  c,  of  the  scattering  surface, 

gHH(e)  =  (6,  -l)[cos  0  +  (6,  -  sin  0)’^]  (Eq.  2) 

the  coordinate  system  is  chosen  such  that  the  k,  is  in  the  same  plane  as  the  radar  antenna 
look  direction  and  the  internal  wave  travel  direction.  Equation  (1)  can  then  be  used  to 
estimate  the  ocean  wave  spectral  density  at  the  Bragg  wavenumber,  Iq,,  using 

Iq,  =  2  ko  sin  0  (Eq.  3) 

Indeed,  the  cross  section  measurements  show  decreases  up  to  6  dB  and  10  dB  at  L-band 
and  X-band  due  to  damping  of  the  wave  spectral  energy  at  the  Bragg  wavelengths  within 
the  slicks.  The  Bragg  wavenumbers  and  wavelengths  based  on  the  radar  configuration 
used  during  the  October  16  morning  measurements  are  given  in  Table  2.  The  RCS  is 
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TABLE  2. 

Bragg  Wavelength  for  L-Band  and  X-Band  Radar  at  0  =  25° 


L-BAND 

X-BAND 

K 

0.314  rad/cm 

1.96  rad/cm 

K 

0.26  rad/cm 

1.7  rad/cm 

K 

23.7  cm 

3.8  cm 

Wb  * 

16  rad/s 

40  rad/s 

fb 

2.5  Hz 

6.4  Hz 

§  ci)b  is  the  intrinsic  angular  frequency  of  the  Bragg  wave. 
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proportional  to  the  ocean  wave  spectral  density  at  Bragg  wavelengths  24  cm  (L-band)  and 
4  cm  (X-band).  The  relative  attenuation  of  24  cm  waves  compared  to  4  cm  waves  can 
be  estimated  by  forming  the  ratio  of  the  L-band  cross  section  crL°(kb)  to  the  X-band  cross 
section  ax°(kb)  as, 

Ak.L)/a“(k,x)  =  I  g^HH(0)  I  ^  S(W,0)/  I  g^^(Q)  I  ^  S(k,x,0)  (Eq.  4) 

It  can  be  shown  that  the  ratio  |  g‘'HH(0)  MM  S^hh(0)  M  *=  1  when  the  imaginary  part 
of  Im  I  fr  I  »  1 .  We  contend  that  this  approximation  applies  to  the  surfactant  free  and 
surfactant  covered  ocean  when  0  =  25°  based  on  e,  «  73  -85i  (L-band)  and  *  48- 
35i  (X-band).  That  is,  the  skin  depth  (electromagnetic  penetration  depth)  is  large 
compared  to  the  surfactant  microlayer  thickness.  Indeed,  the  skin  depth  of  sea  water  is 
about  1  cm  and  0.2  cm  at  L-band  and  X-band,  respectively,  whereas  the  surfactant 
microlayer  thickness  is  considered  less  than  0.1  cm.  Therefore  the  scattering  coefficient 
gHH(0)  should  not  be  substantially  effected  by  this  microlayer  film.  Based  on  this  result, 
we  approximate  equation  (4)  by, 

AKO/AKx)  »  S(Kl,0)/  S(lq,x,0)  (Eq.  5) 

The  straining  of  the  short  Bragg  waves  by  the  internal  waves  can  also  be  inferred 
using  Bragg  scattering  theory.  The  normalized  differential  cross  section  6(r°/o®  is 
proportional  to  the  product  of  the  internal  wave  horizontal-surface  strain  rate  du/dx  and 
the  relaxation  time  r,  [Alpers,  1985],  given  by, 

6a“/a“  =  -(4-1-7)  Tj  du/dx  (Eq.  6) 


where  7  is  the  ratio  between  the  group  and  phase  velocities  (0.5  for  gravity  waves)  of 


the  Bragg  wave.  The  relaxation  time  depends  on  the  radar  look  direction  relative  to  the 
wind  speed  and  direction,  whereas  the  strain  rate  is  dependent  on  the  internal  wave 
propagation  direction  and  associated  surface  current  gradient  also  relative  to  the  radar 
look  direction.  The  relaxation  time  is  approximately  equal  to  the  reciprocal  wind  wave 
growth  rate  when  the  wind  and  waves  are  traveling  in  the  same  direction  [Plant,  1982] 

T,  =  6  *  ®  25  (u./c)-2  /  cj  (Eq.  7) 

where  u.  is  the  friction  velocity  (about  1/25  of  the  wind  speed  at  10  m  height),  c  is  the 
phase  velocity  and  o)  is  the  angular  frequency  of  the  surface  waves.  Equation  (7)  is,  to 
first  order,  valid  for  a  clean  surface  as  demonstrateo  by  Plant  and  a  surfactant  surface 
as  demonstrated  by  Mitsuyasu  and  Honda  [1986].  Alpers  [1985]  has  estimated  t,  =  39 
sec  for  L-band  Bragg  waves  at  a  wind  speed  of  4  m/s  and  internal  wave  strain  rates  in 
coastal  regions  of  du/dx  *  sec  *.  However,  the  relaxation  times  based  on  equation 
(7)  are  much  smaller  at  shorter  wavelengths,  r,  »  1.6  sec  at  X-band  (4  cm)  for  similar 
wind  speeds.  Therefore  the  predicted  changes  in  the  surfactant  free  surface  spectrum  at 
4  cm  wavelengths  is  a  factor  of  ~25  less  than  the  L-band  modulation  (24  cm).  The 
above  parameters  give. 


da^/o^  ~  18  % 

at  L-band 

(Eq.  8) 

So^/o"  -  0.7  % 

at  X-band. 

(Eq.  9) 

Indeed,  the  cross  sectional  modulation  due  to  straining  of  the  short  waves  is  significantly 
smaller  at  X-band  than  L-band  backscatter.  More  importantly,  the  radar  signature  due 
to  straining  has  been  observed  as  a  series  of  bright  and  dark  bands  in  radar  imagery. 
Equations  (1)  and  (5)  provide  a  method  of  extracting  the  damping  of  24  cm  and 
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4  cm  ocean  surface  waves  and  their  differential  damping  within  a  slick  band.  In  contrast, 
equation  (6)  gives  a  method  of  extracting  the  product  of  the  surface  straining  rate  and  the 
relaxation  time.  Furthermore,  equation  (6)  is  consistent  with  the  bright/dark  bands 
observed  in  radar  imagery  due  to  change  in  sign  of  du/dx.  These  expressions,  equations 
(1)  and  (6),  can  be  used  to  distinguish  between  cross  sectional  radar  (reflectivity) 
modulation  due  to  short  wave  damping  and  straining. 

3.3.  OCEAN  SLICK  RESULTS 

We  investigate  the  slick  related  parameters  for  the  period  10:40-1 1:20  on  October 
16  in  more  detail.  Stacked  plots  of  selected  parameters  in  Figure  5  are  displayed  in  a 
format  simitar  to  Figure  4.  The  cross  section  in  Figure  5  (a)  is  given  in  linear  units, 
rather  than  logarithmic  (dB)  units,  0“dB(kb)  =  10  log  A^b)-  The  ratio  of  the  L-band  to 
X-band  cross  section  is  given  in  Figure  5  (b).  The  cross  sectional  ratio  is  directly  related 
to  the  ratio  of  the  spectral  densities  at  the  24  and  4  cm  wavelengths,  S(ki,L,0)/  S(k(,x,0), 
see  equation  (5).  The  relative  fluorescence  F„  originally  displayed  in  Figure  4  (c)  is  also 
given  in  Figure  5. 

These  limited  measurements  of  surface  slicks  show  several  interesting  features. 
The  relative  fluorescence  F,  is  directly  related  to  the  concentration  of  the  chlorophyll  A 
within  the  slick  bands.  The  ratio  of  24  cm  spectral  damping  to  4  cm  damping  is  the 
largest  within  slick  group  "B".  We  investigated  why  the  L-band  cross  section  shows  a 
greater  backscatter  decrease  within  group  "A"  compared  to  group  "B".  The  wind 
direction  changed  from  60°  relative  to  the  radar  look  direction  in  group  "A"  to  30°  in 
group  "B",  whereas,  the  wind  speed  increased  about  1.5  m/s  from  about  2  m/s  to  about 
3.5  m/s  from  group  "A"  to  group  "B".  The  individual  slick  band  orientations  are  not 
known  but  are  considered  constant  relative  to  the  radar  for  the  two  slick  groups.  The  sea 
surface  temperature  measured  at  15  cm  depth  increased  about  0. 1°  C  from  18.8°C  to 
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Figure  5.  Stacked  Plots  of  Measured  Parameters  for  the  40  Minute  Period  10:40  to 
11:20  PDT  Similar  to  Figure  4.  The  Measured  Parameters  are:  (a)  Radar 
Backscatter  Cross  Section  (o")  in  Linear  Units,  (b)  Ratio  of  the  L-Band  to 
X-Band  Cross  Section,  (c)  Relative  Fluorescence  F^. 


18.9°C  between  group  "A"  and  group  "B". 


One  could  argue  that  the  L-band  backscatter  (24  cm  waves)  might  be  more 
sensitive  to  the  1.5  m/s  wind  speed  increase  than  the  X-band  backscatter  (4  cm  waves). 
However,  this  is  unlikely  since  4  cm  waves  respond  to  wind  speed  changes  more  rapidly 
than  24  cm  waves.  A  more  likely  explanation  is  that  the  L-band  backscatter  is  more 
sensitive  to  the  wind  direction  change  than  the  X-band  backscatter.  That  is,  the 
horizontally  polarized  Bragg  scattering  angular  wind  dependence  is  more  sensitive  to  24 
cm  waves  than  4  cm  waves.  We  expect  the  largest  damping  when  the  wind  is  blowing 
along  the  radar  look  direction  for  X-band  (4  cm). 

The  largest  values  of  relative  fluorescence  (F,  ®  7)  near  10:47  and  11:10  are 
approximately  co-located  with  the  minimum  values  in  the  cross  section  for  both  L-band 
(24  cm  waves)  and  X-band  (4  cm  waves).  However,  some  of  the  locations  of  the 
individual  slicks  do  not  exactly  coincide  when  the  cross  section  decreases.  We  did  not 
investigate  these  location  discrepancies  further  in  the  present  work. 

The  slick  characterization  as  a  function  of  changing  wind  speed  and/or  direction 
is  of  interest.  A  considerable  amount  of  radar  imagery  of  surface  slicks  has  been  acquired 
in  recent  years.  However,  these  synthetic  aperture  radar  (SAR)  images  are  instantaneous 
two-dimensional  snapshot  of  the  slicks,  an  example  is  given  in  Figure  6.  That  is,  changes 
in  cross  section  as  a  function  of  wind  speed  or  direction  are  difficult  to  determine  using 
SAR  imagery.  Ship-board  X-band  radar  modulations  have  been  reported  similar  to 
aircraft  X-band  SAR  measurements  when  internal  waves  travel  near  the  range  direction 
of  the  radars;  see,  [Shuchman  et  al.,  1988].  However,  for  other  directions  smaller  image 
modulations  were  observed.  Indeed,  the  lack  of  bright/dark  type  radar  signatures  in  the 
present  measurements  could  possibly  be  due  to  the  orientation  of  the  slicks  relative  to  the 
radar  look  direction  or  alternately,  slicks  not  associated  with  internal  waves.  We  expect 
the  maximum  bright/dark  signature  when  the  slick  groups  are  traveling  in  the  same 
direction  as  the  radar  look  direction.  The  cross  sectional  modulation  due  to  straining  of 
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Aircraft  Flight  Direction 


Figure  6.  Synthetic  Aperture  Radar  Image  of  Internal  Wave  Slicks  in  the  Georgia 
Straits.  This  L-Band  Radar  Image  Shows  Long  Bright  and  Dark  Bands 
Induced  by  a  Group  of  Internal  Waves  Traveling  Toward  the  Top  of  the 
Image.  The  Bright  Return  Near  the  Rear  of  the  Wave  Group  is  a  Ship. 
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the  short  waves  at  L-band  is  in  reasonable  agreement  with  the  Bragg  scattering  theory. 
However,  theory  predicts  that  the  X-band  modulation  is  typically  a  factor  of  ~  25 
smaller,  see  equations  (8)  and  (9).  Observed  radar  image  modulations  show  similar 
modulations  at  both  L-band  and  X-band;  see,  e.  g.,  Lyzenga  and  Bennett  [1988]. 

The  dissipation  of  wind  driven  short  waves  is  significantly  smaller  for  the 
surfactant  free  sea  than  within  a  slick.  The  surface,  instead  of  being  free,  is  to  a  first 
approximation  represented  by  an  inextensible  (incompressible)  microlayer  within  the 
slick.  This  is  equivalent  to  the  horizontal  velocity  of  the  particles  vanishing  at  the 
interface  between  the  microlayer  and  the  sea  water  sub-layer  with  the  vertical  gradient 
of  the  velocity  directly  related  to  the  losses  due  to  viscosity  at  the  interface.  The  wind 
input  to  a  surfactant  free  sea  must  be  sufficient  to  overcome  the  wind  speed  threshold 
which  dependents  on  bulk  viscosity  of  the  sea.  This  threshold  is  about  1  m/s  which 
means  these  short  waves  are  essentially  always  present  on  the  open  (clean)  ocean.  The 
exponential  decay  time  due  to  the  viscosity,  t,  of  the  short  gravity  waves  based  on  a 
surfactant  free  surface  (constant  surface  tension)  is, 

T  *  Xb^/87r^i'  (clean  surface)  (10) 

where  the  spectral  density  of  the  short  waves  decreases  as  exp(-  2\Jt),  v  is  the  kinematic 
viscosity  for  sea  water,  v  =  1.78  x  10’^  cm^  s  '  then  r  »  0.712  where  the  Bragg 
wavelength  is  given  in  centimeters.  We  obtain  for  the  4  cm  and  24  cm  waves,  * 
400  s,  Tf4  as  10  s  which  indicates  that  the  shorter  (4  cm)  gravity  waves  are  more  rapidly 
damped  (shorter  decay  time)  than  the  24  cm  waves.  However,  these  waves  are 
continuously  generated  by  the  wind  input.  Therefore  both  the  4  and  24  cm  waves  are 
always  present  even  for  light  winds. 

The  viscous  damping  within  slicks  is  much  larger.  We  consider  for  the  moment 
that  the  short  waves  as  generated  in  adjacent  surfactant  free  regions  next  to  the  slicks  and 
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propagated  into  the  slick  region  where  they  are  rapidly  attenuated  resulting  in  the  cross 
section  decreases  as  shown  in  Figures  (4),  (5)  and  (6).  The  exponential  decay  time,  Tf, 
for  the  short  waves  within  a  slick  where  large  concentrations  of  organic  material  are 
present  is  approximated  by  an  inextensible  microlayer  as  given  by  Lamb  [1932], 

Tf  *  (S/vcjbkb^)*'^  (inextensible  microlayer)  (11) 

The  decay  time  for  short  gravity  waves  are  then  given  by,  tj  »  0.375  X,,  based  on 
deep  water  dispersion  where  is  given  in  centimeters.  The  shorter  wavelengths  are 
damped  more  strongly  than  the  longer  wavelengths,  (XbL/Xbx)  =  9.86.  However,  the 
wind  input  must  also  be  considered  before  the  radar  extracted  damping  can  be  compared 
with  equation  (11). 

Applying  equation  (11)  to  the  24  cm  and  4  cm  waves,  we  obtain  tr4  »  20  s,  Tf4 
«  2  s.  The  viscous  decay  time  for  both  the  free  surface  and  the  inextensible  microlayer 
as  a  function  of  Bragg  wavelength  is  given  in  Figure  7.  We  have  assumed  for  the 
moment  that  no  short  waves  are  generated  within  the  slick.  Therefore,  the  distance  the 
short  waves  propagate  into  an  inextensible  microlayer  depends  on  the  group  velocity  of 
the  short  gravity  waves,  about  6  m  for  24  cm  waves  and  about  0.2  m  for  4  cm  waves. 
This  idealized  model  is  based  on  a  slick  which  changes  abruptly  from  a  surfactant  free 
region  to  an  inextensible  microlayer  with  wave  reflection  at  the  boundary  neglected. 
However,  the  wind  input  inside  slicks  is  required  before  a  realistic  representation  of  the 
waves  inside  the  slick  bands  can  be  modeled.  Intuitively,  a  slick  band  should  approach 
an  inextensible  microlayer  only  near  its’  center  with  a  transition  to  a  relatively  surfactant 
free  surface  with  distance  away  from  the  center.  The  microlayer  is  essentially  elastic  in 
this  transition  region.  A  realistic  representation  of  the  damping  must  also  take  into 
account  the  concentration.  We  investigate  the  wind  input  and  the  wave  damping  in  terms 
of  the  ERIM  ocean  model  below. 
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Figure  7.  The  Exponential  Decay  Time  for  Both  the  Free  Surface  and  the  Inextensible 
Microlayer  as  a  Function  of  Bragg  Wavelength  Between  4  cm  and  24  cm. 
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3.4  EOM  MODEL 


The  ERIM  Ocean  Model  is  an  integrated  hydrodynamic  and  electromagnetic 
model  resulting  from  the  short  ocean  wave  generation  and  dissipation  within  the  marine 
microlayer  [Tanis  et  al.,  1989].  The  generation  and  dissipation  include  both  wave 
damping  and  straining  of  the  short  waves.  We  use  this  model  to  investigate  the  wave 
damping  since  wave  straining  has  already  been  considered. 

The  damping  of  the  short  waves  has  been  formulated  in  terms  of  the  action 
density  equation  including  wind  forcing  and  viscous  damping.  The  RCS  is  calculated  at 
the  wavelength  corresponding  to  the  Bragg  resonance.  This  modeled  cross  section  given 
in  terms  of  the  natural  surfactant  concentration,  or  more  precisely  the  film  pressure.  The 
one-dimensional  film  pressure  r  is  defined  as  the  force  per  unit  length  exerted  by  the 
film  as  it  spreads  horizontally  across  the  ocean  surface.  The  model  input  range  of  film 
pressure  is  0.01  -  32  dynes/cm.  The  film  pressure  for  a  clean  (surfactant  free)  ocean 
surface  corresponds  to  approximately  0,01  dynes/cm  and  the  concentration  of  the 
surfactant  near  the  center  of  the  slick  band  is  thought  to  be  on  the  order  of  16  -  32 
dynes/cm.  The  variation  of  RCS  as  a  function  of  operating  frequency  (L-,  C-,  and  X- 
band)  and  film  pressure  are  shown  in  Figure  8.  The  polarization  and  incidence  angle  0 
were  selected  to  correspond  to  SLIX  ’89,  namely  vertical  transmit  and  vertical  receive 
and  0  =  25°.  The  resulting  plots  are  not  very  sensitive  to  the  polarization  or  the 
incidence  angle.  The  radar  operating  wavelength  includes  the  dual  wavelength 
configuration  used  during  SLIX  ’89  (  24  cm  and  4  cm  ocean  wavelengths)  and  a  C-band 
(5.3  GHz),  an  intermediate  ocean  wavelength,  6  cm. 

The  EOM  model  results  show  little  change  in  cross  section  at  L-band  (24  cm), 
significant  changes  at  C-band  (6  cm)  and  substantial  change  at  X-band  (4  cm).  The  cross 
section  shows  little  change  with  radar  wavelength  when  x  =  0.01  or  1.0  dyne/cm.  The 
cross  section  was  essentially  constant  with  radar  wavelength  except  for  a  decrease  of 
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about  6  dB  at  X-band  when  t  =  2.0  dynes/cm.  The  cross  section  shows  significant 
attenuation  at  C-band  and  X-band  for  larger  film  pressures,  x  >  4.0  dynes/cm  with  the 
largest  decrease  at  C-band  rather  than  at  the  short  radar  wavelength  X-band.  We  believe 
this  larger  damping  at  C-band  than  L-band  or  X-band  is  due  to  ocean  wave  coupling 
between  the  traverse  and  longitudinal  wave;  see,  e.  g.  Alpers  and  Huhnerfuss  [1989]. 

The  dampening  of  the  short  waves  is  greater  at  4  cm  compared  to  24  cm  as 
measured  across  the  slicks  and  as  calculated  by  the  EOM  model.  However,  the  model 
shows  a  damping  ratio  of  about  7  -  8  dB  whereas  the  measurements  give  a  larger 
damping  ratio  about  12  dB,  a  factor  of  16.  Therefore  the  EOM  model  underestimates  the 
damping  at  4  cm  compared  to  24  cm  by  about  4  dB.  We  suggest  that  if  more 
comprehensive  measurements  are  taken,  the  validity  of  EOM  may  be  determined  which 
might  include  a  cross  section  dependence  on  the  relative  angle  between  the  radar  look 
direction  and  the  wind  and  internal  wave  directions.  A  validated  EOM  model  could  be 
used  to  extract  the  film  pressure  from  the  comparison  of  the  model  and  measurement. 
The  film  pressure  is  a  measure  of  the  concentration  of  the  surfactant  material. 
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4.0.  RADIOMETRIC  MODEL 


The  radiometric  temperature  Tg  for  a  near  surface  down-looking  radiometer  over 
a  rough  half  space  (wind  sea)  based  on  thermodynamic  equilibrium  is  related  to  the  sea 
surface  temperature  Ts  (in  °K)  by, 


(12) 


where  S  is  the  salinity  of  the  ocean,  U  is  the  wind  speed  at  20  m,  e  is  the  emissivity  of 
the  rough  ocean,  and  Ts(°C)  =  Ts(°K)  -  273.16.  The  emissivity  for  a  smooth  sea  at 
normal  incidence  is  given  by. 


=  1 


/ 


/ 


(13) 


where  f  is  the  complex  dielectric  constant  of  sea  water. 

The  dielectric  constant  e  of  sea  water  is  dependent  on  the  microwave  frequency 
0)  in  rad/sec.,  ionic  conductivity  a  in  mhos/m,  relaxation  constant  t  in  sec.,  and  static 
dielectric  constant  c,  [Debye,  1929]  as  given  by, 


e 


es  - 


1  (jo  t) 


1  - 


(14) 


where  =  8.854  x  10  '^  is  the  permittivity  of  free  space  in  farads/m,  is  the 
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permittivity  at  a  high  frequency,  oj  =  2xf  is  the  radian  frequency  with  f  in  Hz,  a  is  an 
empirical  constant  and  a,  t,  and  e  are  all  functions  of  S  and  Tj,.  Empirical  values  of  a, 
T,  and  e  as  a  function  of  S  and  Tj  were  given  by  Klein  and  Swift  [1977]  for  the 
microwave  frequencies  1.43  and  2.65  GHz.  We  use  these  empirical  values  to  investigate 
a  suitable  radiometric  frequency  for  which  the  emissivity  dependence  on  salinity  is 
negligible,  see  section  3.1.  The  selection  of  this  frequency  means  that  equation  (12)  can 
be  simplified  with  Tg  depending  on  only  U  and  Tg, 


=  e(U.T,)Ts. 


(15) 


We  now,  in  principle,  can  extract  the  wind  speed  U  and  surface  temperate  Tj  using 
equation  (15)  and. 


U  =  fUoo) 


(16) 


where  is  the  scatterometer  backscattering  cross  section  and  f  '  is  the  inverse  of  the 
function  f  given  by. 


^0  =  f(U)  =  aU^ 


(17) 


where  a  and  b  are  empirical  constants  dependent  on  the  scatterometer  frequency  and  the 
wind  direction.  The  extracted  winds  at  20  m  are  accurate  to  about  ±  1  m/s  based  on 
aircraft  circle  flights,  upwind,  crosswind  and  downwind;  see,  e.  g.,  Donelan  and  Pierson 
[1987].  If  the  effect  of  boundary  layer  stability  is  removed,  increased  accuracy  is 
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expected. 

The  emissivity  dependence  on  U  and  Tj,  given  in  equation  (15),  has  been  modeled 
quantitatively  by  Wilheit  [1979].  This  numerical  model  originally  developed  by  Stogryn 
[1967]  has  been  simplified  by  Wilheit  based  on  an  isotropic  normal  distribution  for  the 
ocean  wave  slopes.  Wilheit  obtained  good  agreement  between  the  model  and  radiometer 
measurements  based  on  this  simplification  with  the  variance  of  the  slope  increasing 
linearly  with  the  radiometric  frequency.  Furthermore,  for  high  wind  speeds,  he  also 
included  the  effect  of  foam  on  the  radiometric  temperature.  He  treated  foam  as  partially 
obscuring  the  surface  independent  of  polarization.  We  use  this  model  in  the  present 
work.  Equation  (15)  and  this  model  are  used  to  estimate  the  emissivity  dependence  on 
U  and  Ts. 


4.1.  SEA  SURFACE  TEMPERATURE  RESULTS 


The  optimum  radiometric  frequency  is  selected  such  that  the  emissivity 
dependence  on  salinity  is  negligible  and  also  such  that  the  atmospheric  effects  due  to 
water  vapor  and  clouds  are  minimized.  This  results  in  the  lowest  frequency  consistent 
with  negligible  dependence  on  salinity.  Blume,  et  al.  [1978]  used  radiometers  operating 
at  1.43  and  2.65  GHz  to  measure  salinity  with  aircraft  flights  over  an  estuary.  Since  both 
of  these  frequencies  are  sensitive  to  salinity,  an  optimum  frequency  must  be  larger  than 
2.65  GHz.  We  use  equations  (13)  and  (14)  and  Klein  and  Swift’s  results  in  the  range  3 
<  f  <  10  GHz  to  investigate  the  minimum  frequency  which  results  in  an  error  in 
radiometric  temperature  of  less  than  0.04°C  per  ppt  in  salinity.  The  salinity  of  surface 
water  remains  relatively  constant  for  long  periods  of  time  (near  34.6  ppt)  with  a 
geographic  latitude  variation  of  about  1  ppt  provided  isolated  areas  are  excluded.  These 
areas  include  coastal  regions,  especially  estuaries  and  melting  ice  in  polar  regions 
[Williams,  1962]. 
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The  results  of  our  investigation  are  shown  in  Figure  9  (a)  with  the  emissivity  as 
a  function  of  Tg  with  an  optimum  frequency  of  6  GHz;  while,  Figure  9  (b)  illustrates  the 
same  variation  at  2.65  GHz  for  comparison.  The  range  of  salinity  was  varied  from  20 
ppt  to  40  ppt  in  intervals  of  2  ppt.  Each  ppt  interval  change  in  Figure  9  (a)  corresponds 
to  <  0.00013  change  in  emissivity  or  equivalently  <  0.038  °C/ppt.  Each  ppt  interval 
in  Figure  9  (b)  corresponds  to  <  0.0007  or  <  0.22  °C/ppt.  These  results  are  based  on 
Ts  <  20  °C.  Therefore  the  optimum  radiometric  frequency  is  near  6  GHz  as  given  in 
Table  4  (a).  Indeed,  a  typical  salinity  change  of  1  ppt  with  latitude  causes  a  change  in 
Ts  of  0.038  °C  at  C-band  (6  Ghz)  and  0.22  °C  at  S-band  (2.65  GHz)  for  normal 
incidence. 

We  now  consider  the  emissivity  dependence  on  only  U  and  Tg  using  equation  (15) 
and  Wilheit’s  model.  The  variation  of  the  emissivity  with  incidence  angle  for  three  wind 
speeds,  U  =  6,  10  and  20  m/s  are  given  in  Figure  10  (a).  Wind  speeds  less  than  6  m/s 
have  negligible  foam  cover.  The  upper  three  curves  correspond  to  vertical  polarization 
and  the  lower  three  curves  to  horizontal  polarization.  The  vertically  polarized  curves  tend 
to  intersect  at  an  incidence  angle  0  between  50°  and  60°.  Therefore  these  incidence 
angles  correspond  to  a  minimum  dependence  on  wind  speed  which  is  the  same  result  as 
previously  reported  by  Stogryn  [1967].  The  percent  of  effective  foam  cover  for  wind 
speeds  from  8  to  20  m/s  based  on  the  model  are  given  in  Table  3.  Figure  10  (b)  shows 
the  same  results  as  Figure  10  (a)  except  the  differential  emissivity  is  displayed  on  the 
vertical  axis.  This  normalized  differential  emissivity  6e/e,  in  percent  is  defined  as, 


6e(d)  ^  (es(0)  -  e(d)) 

es(d) 


(18) 


where  e,(0)  is  here  the  emissivity  of  the  smooth  sea  as  a  function  of  incidence  angle  and 
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TABLE  3. 

Wind  Speed  and  Associated  Percent  Effective 
Foam  Cover  Based  on  Wilheit’s  Model 


U(m/s) 

FOAM(%) 

6 

0 

8 

.33 

10 

.99 

12 

1.65 

14 

2.31 

16 

2.97 

18 

3.63 

20 

4.31 
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TABLE  4. 

Radiometric  Brightness  Temperature  Sensitivity  to: 

(a)  Salinity  at  Normal  Incidence,  (b)  Wind  Speed,  S  =  34  ppt  And  0  =  60° 


(a) 


OPERATING  FREQUENCY 

C-BAND  (6  GHz) 

S-BAND  (2.65  GHz) 

ATb/AS 

0.038°C/ppt 

0.22°C/ppt 

(b) 


C-BAND  (6  GHz) 

POLARIZATION 

VERTICAL 

HORIZONTAL 

ATb/AU 

<  0.12°C/m/s 

<  2.0°C/m/s 
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e,H(0)  and  e,v(0)  are  the  corresponding  vertical  and  horizontal  polarizations  of  e,(0) 
given  by, 


cos 


cos 


6  -  /e  -  sin^  6 
d  +  \/e  -  sin^d 


2 


(19) 


eje)  =  1 


/ 


e  cos  6  -  \/e  -  sin^d 
e  cosd  +  \/e  -  silled 


I 


2 


(20) 


The  three  lower  curves  correspond  to  vertical  polarization  and  the  three  upper  curves 
correspond  to  horizontal  polarization.  The  differential  emissivity  for  the  vertical 
polarization  consistently  shows  a  lower  percent  variation  as  a  function  of  incidence  angle 
and  wind  speed. 

We  next  consider  the  emissivity  dependence  on  sea  surface  temperature  for  a 
range  of  wind  speeds,  U  =  6-20  m/s  at  0  =  60°.  Figure  11  (a)  shows  the  emissivity 
dependence  on  U  and  Tj  for  horizontal  polarization  whereas  Figure  11  (b)  shows  this 
same  dependence  for  vertical  polarization.  The  primary  result  here  is  that  the  vertical 
polarized  Tb''  is  only  weakly  dependent  on  wind  speed  when  6  m/s  <  U  <  10  m/s 
whereas  at  higher  wind  speeds  the  change  in  emissivity  is  Ae/AU  <  0.00042/m/s  or  in 
terms  of  Tg'^  is  ATbVAU  <  0.12'’C/m/s.  In  contrast,  the  change  in  emissivity  for 
horizontal  polarization  is  Ae/AU  ^  0.007/m/s  or  ATb”/AU  ^  2.0°C/m/s  when  U  >  6 
m/s.  These  Tg  sensitivities  to  wind  speed  are  summarized  in  Table  4  (b).  The  maximum 
error  is  ±  0. 12 °C  and  ±  2.0°C  for  Tg''  and  Tg”  respectively  based  on  the  scatterometer 
wind  accuracy  of  ±  1  m/s.  That  is,  this  analysis  assumes  no  errors  in  the  model  and  the 
measurement  of  the  Tg''  and  Tg”.  This  result  suggests  that  a  scatterometer  can  be  used 
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Figure  1 1 .  Emissivity  Variation  with  Sea  Surface  Temperature  Parameterized  in  Surface 
Wind  SpMd,  6  -  20  m/s.  Incidence  Angle  =  60“  and  Salinity  =  34  ppt;  (a) 
Horizontal  Polarization,  and  (b)  Vertical  Polarization. 
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to  correct  a  radiometer  with  a  vertical  polarized  pencil  antenna  beam  operating  near  6 
GHz  with  an  incidence  angle  near  60°. 
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5.0  PROPOSED  FURTHER  ANALYSIS  AND  FUTURE  EXPERIMENTS 

We  propose  analysis  of  the  radiometer  measurements  on  October  16  during 
SLIX  ’89  when  a  temperature  front  was  intersected  by  the  research  vessel,  see  section 
5.1.  Furthermore,  we  propose  two  new  experimental  programs; 

(1)  Simultaneous  scatterometer,  radiometer  and  SCUMS  measurements  from  a 
ship  platform  across  natural  slicks  at  times  when  radar  images  of  the  slick 
are  available.  This  experiment  is  similar  to  SLIX  ’89  except  we  require  radar 
imagery  in  addition  to  the  other  measurements. 

(2)  Simultaneous  scatterometer  and  radiometer  measurements  from  an  ocean 
tower  acquired  over  a  sufficiently  long  period  of  time  that  a  wide  range  of 
ocean/atmospheric  conditions  are  observed  including  ocean  temperature 
fronts  and  wind  speeds  ranging  over,  say,  0  to  20  m/s. 

5.1  RECC.^/f'-lENDED  ADDITIONAL  ANALYSIS  OF  SLIX  ’89 

The  research  vessel  Wecoma,  during  SLIX  ’89,  traversed  a  temperature  front 
where  the  temperature  change  was  about  1°  C  to  2°  C  on  October  16.  We  suggest  the 
radiometric  temperatures  be  analyzed  for  the  passive  microwave  sensors  on  both  sides 
of  this  temperature  front. 

5.2  PROPOSED  FUTURE  EXPERIMENTS 

Two  experiments  are  proposed  as  summarized  in  section  5.0. 

The  first  experiment  would  include  simultaneous  scatterometer,  radiometer,  and  SCUMS 
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measurements  across  natural  occurring  slicks  when  auxiliary  quantitative  information  is 
available  on  both  the  surface  wind  speed  and  direction,  and  the  slick  spacing  and 
orientation.  Second,  we  propose  an  experiment  to  establish  the  feasibility  of  extracting 
sea  surface  temperature  from  simultaneous  active/passive  measurements. 

The  slick  experiment  would  require  an  imaging  radar  be  flown  over  the 
experiment  area  to  measure  the  slick  spacing  and  orientation  in  addition  to  the  microwave 


and  SCUMS  measurements.  This  imagery  could  also  be  useful  in  establishing  the 
mechanism  responsible  for  the  slicks  as  well  as  any  azimuthal  dependence  of  the  slick 
orientation  relative  to  the  scatterometer  look  direction.  Furthermore,  the  radar  images 
may  be  useful  in  spatially  locating  any  surface  wind  fronts  in  the  area.  The  anemometer 
measured  winds  would  be  continuously  acquired  aboard  the  research  vessel. 

The  proposed  experiment  to  extract  the  sea  surface  temperature  would  use  a 
surface-based  radiometer  operating  at  6  GHz  with  data  acquired  on  both  horizontal  and 
vertical  polarized  antennas  using  a  pencil  beam,  say  A6  <  2°,  selected  at  an  incidence 
angle  near  55°  as  suggested  by  Figure  10  (b).  These  measurements  would  be  acquired 
in  conjunction  with  scatterometer  and  wind  anemometer  measurements  in  an  area  where 
high  winds  and  ocean  temperature  fronts  frequently  occur  allowing  the  validation  of  the 
model  and  the  measuring  technique. 
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6.0  SUMMARY 

The  limited  measurements  reported  here  show  that  the  peaks  in  fluorescence, 
related  to  the  peak  concentration  of  Chlorophyll  A,  are  co-located  with  the  damping  of 
the  short  gravity  waves  at  both  4  cm  and  24  cm  wavelengths.  The  maximum  damping 
ratio  between  the  24  cm  and  4  cm  waves  is  at  least  16  for  a  wind  direction  of  about  30° 
with  respect  to  the  radar  look  direction  and  wind  speeds,  3-4  m/s. 

Imaging-radar  slick  signatures  displaying  narrow  bright  bands  on  the  leading  edge 
of  a  series  of  wider  dark  slick  bands  have  been  interpreted  as  due  to  short  wave  straining 
by  internal  waves.  We  assert  that  the  wide  dark  bands  following  these  bright  bands,  when 
observed,  are  primarily  due  to  wave  damping.  That  is,  these  dark  bands  are  related  to 
increased  concentration  of  the  marine  surface  microlayer  material.  Therefore,  the  radar 
roughness  of  slicks  on  some  occasions  may  be  caused  by  both  short  wave  damping  and 
wave  straining  of  the  short  gravity  waves. 

We  propose  two  future  experiments:  (1)  simultaneous  shipboard  microwave 
sensor  and  SCUMS  measurements  across  naturally  occurring  slicks  similar  to  SLIX  ’89 
except  that  radar  imagery  would  also  be  acquired  during  the  experiment,  and  (2) 
simultaneous  active/passive  measurements  from  an  ocean  tower  using  a  radiometer  with 
both  vertical  and  horizontal  polarization,  pencil  beam  antenna  patterns  at  an  incident 
angle  of  55°  and  an  operating  frequency  of  6  GHz.  In  addition  to  these  radiometer 
measurements,  simultaneous  scatterometer  and  anemometer  measurements  would  also  be 
acquired  on  the  tower. 
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ABSTRACT 

During  the  October  1989  Marine  Hicrolayer 
Slick  Experiment  (SLIX'89)  which  occurred  of'f  the 
coast  of  San  Diego  and  Santa  Rosa  Island,  L-.  C-, 
X-band  scatterometer  and  C-  and  Ka-band  radiometer 
measurements  were  made  from  the  R/V  Wecoma  of 
surface  slicks  associated  with  fronts  and 
naturally  occurring  internal  waves  (IWs).  The 
Wecoma  operated  jji  situ  measuring  devices  as  it 
traversed  the  IW  generated  slicks.  The 
instrumentation  which  operated  coincidently  with 
the  active  and  passive  microwave  sensors  included 
the  SCUMS  catamaran  with  c  jrface  potential,  laser 
slope,  ADCP,  thermistor  f  i,  and  second  harmonic 
surface  probes. 

Microwave  measurement  were  made  both  in  the 
slick  and  surfactant  free  areas  for  comparison 
purposes.  A  range  of  inc  dent  angles  were 
collected  to  demonstrate  llragg  resonant  effects 
due  to  the  slick  material  Initial  examination 
shows  significant  di  fferciice ;  in  both  radar 
backscatter  (cr°)  and  microwave  brightness 
temperature  due  to  the  presence  of  slicks. 

INTRODUCTION 

The  ability  to  discriminate  ocean  surface 
features  caused  by  slicks,  fronts,  and  eddies,  the 
microwave  response,  backscatter  and  emission,  of  a 
wind  driven  sea,  were  thi'  focus  of  this 
investigation.  Overall  goals  include  improving 
the  estimates  of  windspeud  and  direction  and  sea 
surface  temperature  when  using  data  products 
acquired  via  satellite  sensors. 

MEASUREMENTS 

The  microwave  and  slick-ocean  measurements 
(See  Table  1)  were  designed  to  retrieve  inform¬ 
ation  from  patches  of  ocean  surface  observed 
coincidently  with  each  sensing  system.  During  the 
measurement  period  winds  ranged  from  2  to  15 
knots,  with  wave  classes  ranging  from  glass-like 
to  white-capped.  Slicks  occurred  in  distinct 
narrow  bands  or  were  regional.  Fronts  with  1° 
temperature  difference  and  ocean  eddies  were 
observed  during  transit. 

Of  particular  inter(^5t  is  the  relationship 
of  slick  properties  to  the  microwave  response  as 
examined  in  the  frequency  jiid  angular  response 
characteristics  of  the  backi.eatter  depression 
response  determined  from  non-slick  and  slick 
conditions.  SI  ick-induceil  I'oughness  variations 
were  directly  sensed  usimi  the  scatterometers. 
Radiometer  measurements  were  made  to  determine  the 
resultant  effect  on  einisi  vi-y. 


Additionally,  the  hypothesis  that 
backscattering  is  explained  by  Bragg  scattering 
was  tested  by  comparing  radar  frequency-angle 
combinations  which  have  the  same  Bragg  wavelength 
(i.e.  X-20°  8  C-40°,  C-20°  8  L-70°). 

EXAMPLE  RESULTS 

In  the  example  shown  in  Figure  1,  the 
specular  backscatter  component  was  enhanced  at 
angles  near  vertical  under  slick  conditions,  while 
at  angles  from  40°  to  60°,  enhancements  to  5  dB 
were  produced  during  non-slick  conditions  at  X- 
band. 
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Figure  1.  Backscatter  Depression  Versus  Incidence 
Angle  at  X-,  C-,  and  L-band  and  VV- 
Polarization  During  Slick  and  Non-Slick 
Conditions. 

Table  1.  SLIX  '89  Sensor  Description 

Sensors  Measurements 

Radars  Scattering  Coefficients 

1.5,  5,  10  Gllz 
10°  to  85°  Incidence 

Radiometers  Briglitncss  Temperatures 

Emissivities 

5.6  and  35  GHz 
10°  to  85°  Incidence 

Chemical  Surface  Potential 

Composition  8  Biology 

Laser  Wave  Slope  Capillary-Gravity  Waves 

ADCP  Near-Surface  Velocity  8 

Turbulence 

Thermistor  Fin  Near-Surface  Temp. 

Depths  of  5,  10,  15  cm 
Microlayer  Sampler  Absorbance  8  Fluorescence 
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ABSTRACT 


Ship-based  radar  scatterometer  measurements  at  L-band  and  X-band  and 
micrclayer  Chlorophyll  fluorescence  measurements  on  the  southern 
California  continental  shelf  were  acquired  across  several  surface  slicks 
during  October,  1989.  These  measurements  allowed  the  simultaneous 
extraction  of  both  the  short  wave  dcunping  at  Bragg  wavelengths  24  cm  and 
4  cm  and  the  Chlorophyll  A  biological  concentration  with  sufficient 
spatial  resolution  to  observe  features  within  individual  slick  bands. 

The  measurements  show:  (1)  the  largest  peaks  in  fluorescence,  a 
measure  of  the  peak  Chlorophyll  A  concentration,  are  co-located  with  the 
short  wave  damping  at  24  cm  and  the  4  cm  wavelengths,  (2)  the  maximum 
damping  ratio,  between  the  24  cm  and  4  cm  waves  is  at  least  16  for  wind 
speeds,  3-4  m/s,  and  (3)  the  radar  roughness  signatures  at  these  short 
Bragg  gravity  wavelengths  suggest  the  slicks  are  caused  by  damping  rather 
than  straining  of  these  waves.  Additional  shipboard  measurements  are 
needed  to  further  quantify  these  preliminary  results. 
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1.0  INTRODUCTION 


Ship-based  radar  backscatter  ( scatterometer )  and  passive  microwave 
(radiometric)  measurements  from  the  ocean  surface  were  acquired  during  the 
October,  1989,  Marine  Micro-layer  Slick  Experiment  (SLIX'89).  Coincident 
in  situ  measurements  were  also  taken  from  a  catamaran,  mounted  with  the 
Self  Contained  Underway  Microlayer  Sampler  (SCUMS),  towed  from  the  side  of 
the  ship  including  a  rotating  drum  onto  which  the  surface  microlayer 
material  was  attached  [Carlson  et  al.,  1988].  This  microlayer  material  was 
collected  from  the  rotating  drum  and  analyzed  using  UV  absorbance  and 
fluorescence  detectors.  These  ship-based  measurements  provided 
simultaneous  high  spatial  resolution  microwave  and  in  situ  measurements 
with  sufficient  resolution  to  identify  coincident  features  inside 
individual  slicks.  The  area  selected  for  these  measurements  was  the 
southern  California  continental  shelf  were  slicks  are  frequently  observed 
[Ewing,  1950]. 

Radar  scatterometers  measure  the  ocean  surface  roughness  at  Bragg 
wavelengths,  3  cm  -  30  cm,  within  the  radar's  antenna  beamwidth. 
Fluorescence  detectors  measure  the  concentration  of  the  biological 
material,  chlorophyll  A;  whereas,  the  UV  absorbance  detectors  (280  nm) 
measure  the  concentration  of  chemical  material.  The  ability  of  the 
microlayer  to  dampen  short  waves  increases  with  the  concentration  of  the 
microlayer  material.  This  concentration  varies  across  slicks  due  to  the 
horizontal  convergence  of  ocean  currents  at  the  surface  and/or  by 
horizontal  convergence  due  to  wind  stress.  That  is,  the  surface  roughness 
near  and  within  slick  bands  vary  with  ocean  current  and  wind.  Slicks  on 
the  southern  California  and  Baja  continental  shelf  can  be  due  to  internal 
waves  currents,  small  scale  ocean  circulation  and  wind  stress.  The 
internal  waves  typically  travel  along  the  thermocline  at  the  emne  speed  as 
the  companion  slick  bands.  These  slicks  are  characterized  by  a  series  of 
slick  packets;  each  packet  separated  by  a  distance  of  15  -  20  km  with  an 
envelope  of  quasi-period  parallel  slick  bands  with  wavelengths  of  about 
300  km  [Ewing,  1950;  and  Apel  and  Gonzalez,  1983].  The  slick 
characteristics  can  be  used  to  distinguish  internal  wave  generated  slicks 
from  other  slick  generating  mechanisms. 

The  surface  roughness  at  short  gravity  wavelengths  associated  with 
quasi-periodic  slicks  have  been  observed  by  a  number  of  workers  using 
different  measuring  methods;  see,  e.  g.,  Ewing  [1950],  LaFond  [1962], 
Rufenach  and  Smith  [1985],  Alpers  [1985],  Shuchman  et  al.  [1988],  Lyzenga 
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and  Bennett  [1988]  and  Gasparovic  et  al.  [1988].  Alternating  bright  and 
dark  bands  in  radar  images  have  been  interpreted  as  horizontal  current 
gradients  produced  by  internal  waves;  see,  e.  g.,  Alpers  [1985];  or 
equivalently,  increases  and  decreases  in  radar  cross-section.  In 
contrast,  radar  cross-section  decreases  are  directly  related  to  the 
damping  of  these  short  waves  due  to  increases  in  the  concentration  of  the 
surface  marine  microlayer  material.  Surface  roughness  changes  have  usually 
been  interpreted  in  terms  of  either  dzunping  or  straining  of  the  short 
waves  but  not  both.  We  believe  that  both  mechanisms  may  be  important  on 
some  occasions. 

The  concentration  of  the  surface  microlayer  material  has  been 
related  to  the  bulk  viscosity  [Carlson,  1982;  Carlson,  1987].  Microlayer 
viscous  deunping  is  caused  by  either  the  bulk  viscosity,  syrupy  nature  of 
the  surface  layer,  or  the  elasticity,  visco-elastic  nature  of  this  thin 
surface  layer.  The  elastic  properties  are  equivalent  to  an  additional 
viscosity  of  the  surface  microlayer  [Levich,  1962].  We  do  not  attempt  to 
distinguish  between  these  two  viscous  damping  mechanisms. 

Detection  of  bands  of  slicks  associated  with  internal  waves  by 
passive  microwave  sensors  have  not  been  reported.  Passive  microwave 
sensors  have  been  primarily  operated  from  spaceborne  platfoms  with  a 
spatial  resolution  >  20  km.  This  resolution  could  explain  the  lack  of 
detection  since  the  slick  packet  spacing  is  typically  10  -  15  km. 
Increases  in  emissivity  have  been  measured  across  manmade  oil  slicks  by 
microwave  sensors  [Bollinger  and  Mannella,  1973]  whereas  one  example  of  a 
large  emissitivity  decrease  has  been  reported  across  a  manmade 
monomolecular  slick  [Hiihnerfuss  and  Alpers,  1983].  Additional  shipboard 
multi-frequency  active  and  passive  microwave  measurements  in  conjunction 
with  in  situ  measurements  across  manmade  and  natural  slicks  are  needed  to 
better  understand  the  radio  perturbations  caused  by  surface  slicks 

This  study  is  limited  to  a  comparison  of  scatterometer  and 
radiometer  measurements  with  in  situ  fluorescence  measurements  near  two 
surface  slick  groups  off  the  California  coast.  These  measurements  are 
interpreted  in  terms  of  straining  and  damping  of  the  short  gravity  waves 
within  these  two  slick  groups. 
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2.0  Experimental  Description 


The  ship-board  SLIX'89  experiment  was  designed  to  traverse  surface 
slicks  off  the  coast  of  southern  California  near  San  Diego  during  a  two 
week  period  in  October  1989.  The  ship  heading  during  the  experiment  was 
selected  to  intersect  slick  groups  visually  observed  from  the  bridge  of 
the  ship.  This  experiment  included  a  variety  of  in  situ  and  remote  sensing 
instruments  including  the  microwave  sensors  and  SCUMS  catamaran  in  situ 
detectors. 

The  microwave  sensors  were  mounted  about  10  m  above  the  sea  surface 
toward  the  bow  of  the  ship,  R/V  Wecoma,  whereas  the  in  situ  instruments 
were  towed  by  a  cateunaran  attached  to  a  boom  on  the  port  side  of  the  ship. 
The  wind  velocity  at  a  height  of  about  10  m  was  measured  near  the  bow  of 
the  ship  and  sea  surface  temperature  was  measured  at  a  depth  of  15  cm.  He 
operated  the  following  microwave  sensors: 

(1)  a  dual  wavelength  radar  scatterometer,  and 

(2)  a  dual  wavelength  radiometer, 

both  mounted  on  the  port  side  with  their  antennas  pointing  perpendicular 
to  the  ship's  track.  The  radar  operating  wavelengths  were  L-band  and 
X-band  (20  cm  and  3.2  cm)  while  the  radiometers  operated  at  C-band  and  Ka- 
band  wavelengths  (5.5  cm  and  0.86  cm).  The  radiometer  outputs  were 
carefully  monitored  for  interference  from  the  scatterometer.  The  relevant 
instrument  parameters  are  given  in  Table  1.  The  radiometers  were 
calibrated  using  a  hot  absorber  and  cold  (near  zenith  sky  scan)  load. 

The  in  situ  measurements  of  interest  were  taken  on  board  the 
catamaran  towed  about  20  m  aft  of  the  microwave  sensors  from  the  port  side 
of  the  ship,  SCUMS  [Carlson  et  al.,  1988].  The  UV  absorbance  is  related  to 
the  concentration  of  the  absorbing  organic  molecules,  whereas  the 
chlorophyll  A  fluorescence  is  related  to  the  concentration  of  the  biota. 
The  SCUMS  measurements  include:  the  fluorescence  and  the  UV  absorption 
measured  in  both  the  microlayer  and  sub-layer.  The  microlayer  material  was 
attached  to  a  rotating  drum  then  skinoned  from  the  drum  surface  whereas  the 
sub-surface  material  was  collected  at  a  depth  of  about  15  cm.  The 
microlayer  is  thought  to  be  50-100  pm  thick.  The  fluorescence  and  UV 
absorbance  show  similiar  variation  across  the  slicks  implying  that  the 
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TABLE  1.  SCATTEROMETER  AND  RADIOMETER  SYSTEM  PARAMETERS 


RJ\ 

lDAR 

RADIO 

METER 

FREQUENCY 

1.5  GHz 

9.38  GHz 

5.4  GHz 

35  GHZ 

WAVELENGTH 

20  cm 

3.2  cm 

5.5  cm 

0.86  cm 

ANTENNA  SPOT  SIZE  at 
e  =  55" 
at  0  =  25"  S 

3.7  m 

1.5  m 

1.3  m 

0.5  m 

5.3  m 
2.1  m 

5.3  m 
2.1m 

INTEGRATION  TIME 

50  sec 

50  sec 

50  sec 

50  sec 

BANDWIDTH 

350  MHz 

575  MHz 

70  MHz 

500  MHz 

POLARIZATION 

II^SBiSESili 

H  or  V 

H  or  V 

RESOLUTION 

1.1  dB 

2.6  dB 

0.21"K 

0.12"K 

ABSOLUTE  ACCURACY 

±  1  dB 

±  1  dB 

±  1"K 

BBS 

*  H  signifies  horizontal  polarization  and  V  signifies  vertical 
polarization. 

$  incidence  angle  used  during  measurements  across  slick  groups  "A"  and  "B” 
on  October  16. 
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concentration  of  the  surface  microlayer  material  can  be  represented  to 
first  order  by  either  detector. 

A  side  view  of  the  ship,  catamaran  and  remote  sensing  instruments 
are  shown  in  Figure  1(a)  and  a  side  view  of  the  catamaran  and  associated 
instrumentation  is  shown  in  Figure  1(b).  A  photograph  of  the  microwave 
antennas  are  given  in  Figure  2.  Selected  measurements  during  the  SLIX 
experiment  are  given  in  the  following  section. 
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Figure  2.  Photograph  of  the  Radio  Microwave  Sensors!  Radar  Scatterometers 
and  Passive  Radiometers. 


3.0  MEASUREMENTS 


Two  surface  slick  groups  off  the  coast  of  San  Diego  during  October 
16  were  selected  for  analysis.  These  two  groups  each  exhibit  three  slick 
bands  during  a  period  when  both  the  microwave  and  in  situ  instruments  were 
operating.  Furthermore,  these  two  slick  groups  were  visually  observed  from 
the  ship  deck  as  recorded  in  the  operational  log.  We  distinguish  between 
groups  and  packets  in  the  present  work:  packets  are  associated  with 
internal  waves;  typically  separated  by  15  -  20  km.  Groups  are  more  closely 
spaced  clusters  of  slick  bands,  several  occurring  within  a  packet,  due  to 
the  complex  nature  of  the  bands  or  generating  mechanisms  other  than 
internal  waves.  The  orientation  of  these  slick  bands  relative  to  the  radar 
look  direction  is  not  known  since  radar  imagery  was  not  available  and  no 
orientation  information  was  recorded  in  the  log.  However,  the  experiment 
objective  was  to  traverse  the  slick  bands.  That  is,  the  fluorescence  and 
radar  measurements  show  a  fluctuating  component  along  the  ship  track. 

We  have  no  direct  measurements  in  the  water  column  that  these  slicks 
were  produced  by  internal  waves  or  any  other  mechanism.  However,  the 
separation  between  the  two  groups  30  min.)  and  a  slick  band  spacing  (3 
-  5  min.),  see  Figure  4,  suggest  spatial  separations;  the  order  of  a  few 
km  and  hundreds  of  meters,  respectively.  The  3-5  min.  (}uasi-per iodic 
spacing  are  consistent  with  internal  wave  produced  slick  bands  provided 
the  ship  speed  (0.5  -  1  m/s)  is  much  larger  than  that  component  of  the 
slick  velocity  along  the  ship  track.  The  spatial  separation  of  a  few  km 
imply  that  the  two  slick  groups  are  not  internal  wave  associated  packets. 
The  grey  areas  indicated  in  Figure  3  give  the  location  of  the  two  slick 
groups  relative  to  the  ship  track  during  October  16.  The  two  slick  groups 
occurred  in  the  morning  hours  between  10:40  to  11:20  Pacific  Daylight  Time 
(PDT) .  There  were  no  clouds  or  rain  during  this  period. 

An  overview  of  the  microwave  and  in  situ  measurements  using  stacked 
plots  of  the  various  parameters  is  given  in  Figure  4.  The  measured 
parameters  (1  min  averages)  are  displayed  on  the  ordinate  and  the  time  in 
minutes  is  displayed  on  the  abscissa  in  all  the  plots.  The  values  of  radar 
cross-section  (c°),  fluorescence  relative  to  the  sub-layer  (F,)  and 
radiometric  brightness  temperature  (Tg)  were  selected  for  display.  The  two 
a°’a  at  L-band  and  X-band  are  indicated  in  Figure  4(a).  The  relative  slick 
fluorescence  (F,)  is  obtained  by  dividing  the  microlayer  fluorescence  by 
the  sub-layer  fluorescence.  The  inherent  delay  in  skimming  and  pumping  the 
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Pacific  Daylight  Time  (PDT) 

91-20067  R1 

Figure  4.  Stacked  Plots  of  Measured  Paranieters  for  the  40  Minute  Period 
10:40  to  11:20  PDT  for  Two  Slick  Wave  Groups  from  R/V  Wecoma, 
Morning  of  October  16,  1989.  The  Measured  Parameters  are:  (a) 
Radar  Backscatter  Cross-Section  [a°) ,  (b)  Slick  Fluorescence 
Relative  to  the  Sub-Layer  (F,)  and  (c)  Radiometric  Brightness 
Temperature  (Tb).  The  Polarization  was  Vertical  for  the  Radiometer 
and  Horizontal  Transmit  and  Horizontal  Receive  for  the  Radar  Both 
at  an  Incidence  Angle  of  6  s  25”. 
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surfactant  microlayer  has  been  corrected.  The  delay  due  to  the  spatial 
separation  of  the  microwave  instruments  and  the  catamaran,  <  20  sec,  has 
not  been  corrected.  The  radiometric  brightness  temperatures  at  C-band  and 
Ka-band  show  no  discernible  change  across  the  slick  bands.  That  is,  any 
surface  slope  decrease  due  to  the  slicks  was  not  detected  by  these 
radiometers. 

The  UV  absorbance  and  the  fluorescence  show  similar  variations 
across  the  two  surface  slick  groups  of  October  16.  We  select  the 
fluorescence  for  detailed  analysis  in  the  present  work.  The  fluorescence 
changes  from  10:40  to  11:20  is  a  measure  of  the  concentration  of  the 
chlorophyll  A  within  the  two  slick  groups  separated  by  a  relatively 
surfactant  free  region.  The  first  group  is  designated  as  group  "A"  and  the 
second  group  is  designated  as  group  "B".  The  three  peaks  in  fluorescence 
F,  within  each  slick  group  are  interpreted  as  three  slick  bands. 

The  relative  fluorescence  F,  increased  to  a  peak  value  of  7  within 
the  slick  group  "A”  near  10:47.  The  radar  cross-section  also  decreasea 
about  10  dB  at  X-band  and  6  dB  at  L-band,  both  near  10:47,  co-located  with 
the  fluorescence  peak.  This  co-location  is  direct  evidence  that  the  short 
Bragg  waves  are  dam.ped  at  the  peak  surfactant  concentration  near  the 
center  of  the  slick  band. 

The  relative  fluorescence  F,  increased  to  a  peak  value  of  8  within 
slick  group  "B"  near  11:10.  Different  radar  cross-section  decreases  are 
apparent  within  slick  group  "B"  compared  to  slick  group  "A".  The  radar 
cross-section  decreased  about  10  dB  at  X-band  similiar  to  slick  group  "A" 
whereas  the  cross  section  decreased  about  2  dB  at  L-band  in  contrast  to 
the  6  dB  decrease  in  slick  group  "A".  Furthermore,  the  X-band  decrease 
occurred  over  a  longer  interval,  about  8  min,  than  group  "A".  This  large 
damping  at  X-band,  reducing  the  backscattered  energy  to  near  the  noise 
floor  of  the  scatterometer,  does  not  allow  co-location  of  the  peak 
fluorescence  F,  with  the  radar  minima  signatures.  That  is,  this  X-band 
damping  occurred  across  several  slick  bands  within  Group  "S",  whereas,  the 
L-band  radar  minima  are  approximately  co-located  with  all  three 
fluorescence  peaks. 

The  wind  speed  changed  from  about  1.5  m/s  to  zd}Out  4  m/s  while  the 
wind  direction  changed  from  about  60"  to  about  30°  with  respect  to  the 
radar  look  direction  from  group  "A"  to  group  "B".  The  wind  was  blowing 
generally  toward  the  south  and  south-east  toward  the  radar  look  direction. 
The  ship  speed  was  about  0.5  m/s.  The  wind  speeds  were  light,  1.5-3  m/s 
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inside  group  "A"  at  a  wind  direction  of  about  60®.  Slightly  higher  wind 
speeds,  3-4  m/s  were  observed  inside  group  "B"  at  a  direction  of  about 
30°.  The  differences  between  the  X-band  radar  cross  section  decreases 
within  the  two  slick  groups  could  be  due  to  the  wind  direction  changes 
between  the  two  slick  groups.  These  differences  are  discussed  in  more 
detail  in  the  following  sections.  Decreases  in  radar  cross  section  co¬ 
located  with  fluorescence  peaks  were  also  observed  within  several  other 
slicks. 
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4.0  BRAGG  SCATTERING  MODEL 


Typical  ocean  remote  sensing  radars  operate  at  microwave  wavelengths 
in  the  centimeter  or  decimeter  range  at  incidence  angles,  9  from  20-70° 
within  which  the  radar  backscatter  is  dominated  by  Bragg  scattering.  The 
radar  cross-section  to  first  order  for  horizontal  transmit  and  receive 
polarization  is  given  by  [Valenzuela,  1978], 

a“(kt)  =  4  ko^  cos^  0  igHH(0)  T  S(k^,0)  (1) 

where  ko  is  the  radar  wavenumber,  S(k,,ky)  is  the  spectral  density  of  the 
ocean  surface  waves  in  rectangular  coordinates,  and  gHH(®)  is  the  first 
order  scattering  coefficient  dependent  on  the  relative  dielectric  constant 
t,  of  the  scattering  surface, 

gHH(0)  =  («,  -l)[cos  0  +  (€,  -  sin  O)'"*]  (2) 

the  coordinate  system  is  chosen  such  that  the  k,  is  in  the  same  plane  as 
the  radar  antenna  look  direction  and  the  internal  wave  travel  direction. 
Equation  (1)  can  then  be  used  to  estimate  the  ocean  wave  spectral  density 
at  the  Bragg  wavenumber,  k*,,  using 

Jq,  =  2  ko  sin  0  (3) 

Indeed,  the  cross-section  measurements  show  decreases  up  to  6  dB  and  10  dB 
at  L-band  and  X-band  due  to  damping  of  the  wave  spectral  energy  at  the 
Bragg  wavelengths  within  the  slicks.  The  Bragg  wavenumbers  and  wavelengths 
based  on  the  radar  configuration  used  during  the  October  16  morning 
measurements  are  given  in  Table  2.  The  radar  cross-section  is  proportional 
to  the  ocean  wave  spectral  density  at  Bragg  wavelengths  24  cm  (L-band)  and 
4  cm  (X-band).  The  relative  attenuation  of  24  cm  waves  compared  to  4  cm 
waves  can  be  estimated  by  forming  the  ratio  of  the  L-band  cross-section 
to  the  X-band  cross-section 

‘7”(kM.)/a'’(k^)  =  |g‘'m(e)  I'  S(k«.,0)/|g’‘m(e)  r  S(k«,0)  (4) 

It  can  be  shown  that  the  ratio  |g‘'HH(®)|‘  /|g’'HH(6)|^  =  1  when  the  imaginary 
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TABLE  2.  BRAGG  WAVELENGTH  FOR  L-BAND  AND  X-BAND  RADAR  AT  6  =  25 


L-BAND 

X-BAND 

K 

0.314  rad/cm 

1.96  rad/cm 

K 

0.26  rad/cm 

1.7  rad/cm 

23.7  cm 

3.8  cm 

^  ‘ 

16  rad/s 

40  rad/s 

2.5  Hz 

6.4  Hz 

S  is  the  intrinsic  angular  frequency  of  the  Bragg  wave. 
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part  of  I(n|e,|  »  1.  We  contend  that  this  approximation  applies  to  the 
surfactant  free  and  surfactant  covered  ocean  when  6  =  25°  based  on  e,  »  73 
-85i  (L-band)  and  e,  ~  48-  35i  (X-band) .  That  is,  the  skin  depth 
(electromagnetic  penetration  depth)  is  large  compared  to  the  surfactant 
microlayer  thickness.  Indeed,  the  skin  depth  of  sea  water  is  eibout  1  cm 
and  0.2  cm  at  L-band  and  X-band,  respectively,  whereas  the  surfactant 
microlayer  thickness  is  considered  less  than  0.1  cm.  Therefore  the 
scattering  coefficient  gHH(^)  should  not  be  substantially  effected  by  this 
microlayer  film.  Based  on  this  result,  we  approximate  equation  (4)  by. 


«^®(k«.)/ff“(k«)  =  S(kfcL,0)/  S(k«,0)  (5) 

The  straining  of  the  short  Bragg  waves  by  the  internal  waves  can 
also  be  inferred  using  Bragg  scattering  theory.  The  normalized 
differential  cross-section  is  proportional  to  the  product  of  the 
internal  wave  horizontal-surface  strain  rate  du/dx  and  the  relaxation  time 
r,  [Alpers,  1985],  given  by, 

=  -(4+7)  r,  du/dx  (6) 

where  7  is  the  ratio  between  the  group  and  phase  velocities  (0.5  for 
gravity  waves)  of  the  Bragg  wave.  The  relaxation  time  depends  on  the  radar 
look  direction  relative  to  the  wind  speed  and  direction,  whereas  the 
strain  rate  is  dependent  on  the  internal  wave  propagation  direction  and 
associated  surface  current  gradienc  also  relative  to  the  radar  look 
direction.  The  relaxation  time  is  approximately  equal  to  the  reciprocal 
wind  wave  growth  rate  when  the  wind  and  waves  are  traveling  in  the  same 
direction  [Plant,  1982]  , 

r,  =  fl'  =  25  (u./c)'^  /  w  (7) 

where  u.  is  the  friction  velocity  (about  1/25  of  the  wind  speed  at  10  m 
height),  c  is  the  phase  velocity  and  u>  is  the  angular  frequency  of  the 
surface  waves.  Equation  (7)  is  ,to  first  order,  valid  for  a  clean  surface 
as  demonstrated  by  Plant  and  a  surfactant  surface  as  demonstrated  by 
Mitsuyasu  and  Honda  [1986].  Alpers  [1985]  has  estimated  r,  =  39  sec  for  L- 
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band  Bragg  waves  at  a  wind  speed  of  4  m/s  and  internal  wave  strain  rates 
in  coastal  regions  of  du/dx  ~  10'^  sec'.  However,  the  relaxation  times  based 
on  equation  (7)  are  much  smaller  at  shorter  wavelengths,  r,  »  1.6  sec  at 
X-band  (4  cm)  for  similiar  wind  speeds.  Therefore  the  predicted  changes  in 
the  surfactant  free  surface  spectrum  at  4  cm  wavelengths  is  a  factor  of 
-25  less  than  the  L-band  modulation  (24  cm).  The  above  parameters  give, 

6c7“/a“  -  18  %  at  L-band  (8) 

6£7®/a®  -  0.7  %  at  X-band.  (9) 

Indeed,  the  cross-sectional  modulation  due  to  straining  of  the  short  waves 
is  significantly  smaller  at  X-band  than  L-band  backscatter.  More 
importantly,  the  radar  signature  due  to  straining  has  been  observed  as  a 
series  of  bright  and  dark  bands  in  radar  imagery. 

Equations  (1)  and  (5)  provide  a  method  of  extracting  the  damping  of 
24  cm  and  4  cm  ocean  surface  waves  and  their  differential  damping  within 
a  slick  band.  In  contrast,  equation  (6)  gives  a  method  of  extracting  the 
product  of  the  surface  straining  rate  and  the  relaxation  time. 
Furthermore,  equation  (6)  is  consistent  with  the  bright/dark  bands 
observed  in  radar  imagery  due  to  change  in  sign  of  du/dx.  These 
expressions,  equations  (1)  and  (6),  can  be  used  to  distinguish  between 
cross-sectional  radar  modulation  due  to  short  wave  damping  and  straining. 
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5.0  DISCUSSION 


We  investigate  the  slick  related  parameters  for  the  period  10:40- 
11:20  on  October  16  in  more  detail.  Stacked  plots  of  selected  parameters 
in  Figure  5  are  displayed  in  a  format  similiar  to  Figure  4.  The  cross- 
section  in  Figure  5  (a)  is  given  in  linear  units,  rather  than  logarithmic 
(dB)  units,  =  10  I09  •  The  ratio  of  the  L-band  to  X-band 

cross-section  is  given  in  Figures  5  (b).  The  cross-sectional  ratio  is 
directly  related  to  the  ratio  of  the  spectral  densities  at  the  24  and  4  cm 
wavelengths,  S(ku.,0)/  S(ku(,0),  see  equation  (5).  The  relative  fluorescence 
F,,  originally  displayed  in  Figure  4  (c)  is  also  given  in  Figure  5. 

These  limited  measurements  of  surface  slicks  show  several 
interesting  features.  The  relative  fluorescence  F,  is  directly  related  to 
the  concentration  of  the  chlorophyll  A  within  the  slick  bands.  The  ratio 
of  24  cm  spectral  damping  to  4  cm  damping  is  the  largest  within  slick 
group  "B”.  We  investigate  why  the  L-band  cross-section  shows  a  greater 
backscatter  decrease  within  group  "A”  compared  to  group  "B”.  The  wind 
direction  changed  from  60°  relative  to  the  radar  look  direction  in  group 
"A”  to  30°  in  group  "B”,  whereas,  the  wind  speed  increased  about  1.5  m/s 
from  about  2  m/s  to  about  3.5  m/s  from  group  "A"  to  group  "B".  The 
individual  slick  band  orientations  are  considered  constant  relative  to  the 
radar  within  the  two  slick  groups.  The  sea  surface  temperature  measured  at 
15  cm  depth  increased  about  0.1°  C  from  18.8°C  to  18.9°C  between  group  "A" 
and  group  "B". 

One  could  argue  that  the  L-band  backscatter  (24  cm  waves)  might  be 
more  sensitive  to  the  1.5  m/s  wind  speed  increase  than  the  X-band 
backscatter  (4  cm  waves).  However,  this  is  unlikely  since  4  cm  waves 
respond  to  wind  speed  changes  more  rapidly  than  24  cm  waves.  A  more  likely 
explanation  is  that  the  L-band  backscatter  is  more  sensitive  to  the  wind 
direction  change  than  the  X-band  backscatter.  That  is,  the  horizontally 
polarized  Bragg  scattering  angular  wind  dependence  is  more  sensitive  to  24 
cm  waves  than  4  cm  waves.  We  expect  the  largest  damping  when  the  wind  in 
blowing  along  the  radar  look  direction  for  X-band 
(4  cm ) . 

The  largest  values  of  relative  fluorescence  (F,  »  7)  near  10:47  and 
11:10  are  approximately  co-located  with  the  minimum  values  in  the  cross- 
section  for  both  L-band  (24  cm  waves)  and  X-band  (4  cm  waves).  However, 
some  of  the  locations  of  the  individual  slicks  do  not  exactly  coincide 
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Figure  5. 
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stacked  Plots  of  Measured  Parameters  for  the  40  Minute  Period 
10:40  to  11:20  PDT  for  Similar  to  Figure  4.  The  Measured 
Parameters  are:  (a)  Radar  Backscatter  Cross-Section  (a")  in  Linear 
Units,  (b)  Ratio  of  the  L-Band  to  X-Band  Cross-Section,  (c) 
Relative  Fluorescence  F,. 
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with  the  cross-section  decreases.  We  do  not  investigate  these  location 
discrepancies  further  in  the  present  work. 

The  slick  characterization  as  a  function  of  changing  wind  speed 
and/or  direction  is  of  interest.  A  considerable  amount  of  radar  imagery  of 
surface  slicks  has  been  acquired  in  recent  years.  However,  these  synthetic 
aperture  radar  (SAR)  images  are  instantaneous  two-dimensional  snapshot  of 
the  slicks,  an  exeunple  is  given  in  Figure  6.  That  is,  changes  in  cross- 
section  as  a  function  of  wind  speed  or  direction  are  difficult  to 
determine  using  SAR  imagery.  Ship-board  X-band  radar  modulations  have  been 
reported  similiar  to  aircraft  X-band  SAR  measurements  when  internal  waves 
travel  near  the  range  direction  of  the  radars;  see,  [Shuchman  et  al., 
1988].  However;  for  other  directions  smaller  image  modulations  were 
observed.  Indeed,  the  lack  of  bright/dark  type  radar  signatures  in  the 
present  measurements  could  possibly  be  due  to  the  orientation  of  the 
slicks  relative  to  the  radar  look  direction  or  alternately,  slicks  not 
associated  with  internal  waves.  We  expect  the  maximum  bright/dark 
signature  when  the  slick  groups  are  traveling  in  the  same  direction  as  the 
radar  look  direction.  The  cross-sectional  modulation  due  to  straining  of 
the  short  waves  at  L-band  is  in  reasonable  agreement  with  the  Bragg 
scattering  theory.  However,  theory  predicts,  that  the  X-band  modulation  is 
typically  a  factor  of  -  25  smaller  ,  see  equations  (8)  and  (9).  Observed 
radar  image  modulations  show  similiar  modulations  at  both  L-band  and  X- 
band;  see,  e.  g.,  Lyzenga  and  Bennett  [1988]. 

The  dissipation  of  wind  driven  short  waves  is  significantly  smaller 
for  the  surfactant  free  sea  than  within  a  slick.  The  surface,  instead  of 
being  free,  is  to  a  first  approximation  represented  by  an  inextensible 
(incompressible)  microlayer  within  the  slick.  This  is  equivalent  to  the 
horizontal  velocity  of  the  particles  vanishing  at  the  interface  between 
the  microlayer  and  the  sea  water  sub-layer  with  the  vertical  gradient  of 
the  velocity  directly  related  to  the  losses  due  to  viscosity  at  the 
interface.  The  wind  input  to  a  surfactant  free  sea  must  be  sufficient  to 
overcome  the  wind  speed  threshold  which  dependents  on  bulk  viscosity  of 
the  sea.  This  threshold  is  about  1  m/s  which  means  these  short  waves  are 
essentially  always  present  on  the  open  ocean.  The  exponential  decay  time 
due  to  the  viscosity,  r,  of  the  short  gravity  waves  based  on  a  surfactant 
free  surface  (constant  surface  tension)  is, 

r  a  X^/Bk^v  (clean  surface)  (10) 
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Aircraft  Flight  Direction  - ► 

Figure  6.  Synthetic  Aperture  Radar  Image  of  Internal  Wave  Slicks  in  the 
Georgia  Straits.  This  L-Band  Radar  Image  Shows  Long  Bright  and 
Dark  Bands  Induced  by  a  Group  of  Internal  Waves  Traveling  Toward 
the  Top  of  the  Image.  The  Bright  Return  Near  the  Rear  of  the  Wave 
Group  is  a  Ship. 
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where  the  spectral  density  of  the  short  waves  decreases  as  exp(-  2t/r),  y 
is  the  kinematic  viscosity  for  sea  water,  v  =  1.78  x  10'^  cm^  s'  then  t  = 
0.712  where  the  Bragg  wavelength  A,,  is  given  in  centimeters.  We  obtain 
for  the  4  cm  and  24  cm  waves,  Tq,  ~  400  s,  r,,  a;  lo  s  which  indicates  that 
the  shorter  (4  cm)  gravity  waves  are  more  rapidly  damped  (shorter  decay 
time)  than  the  24  cm  waves.  However,  these  waves  are  continuously 
generated  by  the  wind  input.  Therefore  both  the  4  and  24  cm  waves  are 
always  present  even  for  light  winds. 

The  viscous  damping  within  slicks  is  much  larger.  We  consider  for 
the  moment  that  the  short  waves  as  generated  in  adjacent  surfactant  free 
regions  next  to  the  slicks  and  propagated  into  the  slick  region  where  they 
are  rapidly  attenuated  resulting  in  the  cross-section  decreases  as  shown 
in  Figures  (4),  (5)  and  (6).  The  exponential  decay  time,  r„  for  the  short 
waves  within  a  slick  where  large  concentrations  of  organic  material  are 
present  is  approximated  by  an  inextensible  microlayer  as  given  by  Lamb 
[1932], 


Tf  «  (8/i/Wt,kb^) (inextensible  microlayer)  (11) 

The  decay  time  for  short  gravity  waves  are  then  given  by,  r,  »  0.375 
based  on  deep  water  dispersion  where  A),  is  given  in  centimeters.  The 
shorter  wavelengths  are  daunped  more  strongly  than  the  longer  wavelengths, 
=  9.86.  However,  the  wind  input  must  also  be  considered  before 
the  radar  extracted  damping  can  be  compared  with  equation  (11). 

Applying  equation  (11)  to  the  24  cm  and  4  cm  waves,  we  obtain, 

~  20  e,  Tf,  s  2  s.  The  viscous  decay  time  for  both  the  free  surface  and  the 
inextensible  microlayer  as  a  function  of  Bragg  wavelength  is  given  in 
Figure  7 .  We  have  assumed  for  the  moment  that  no  short  waves  are  generated 
within  the  slick.  Therefore,  the  distance  the  short  waves  propagate  into 
an  inextensible  microlayer  depends  on  the  group  velocity  of  the  short 
gravity  waves,  about  6  m  for  24  cm  waves  and  about  0.2  m  for  4  cm  waves. 
This  idealized  model  is  based  on  a  slick  which  changes  abruptly  from  a 
surfactant  free  region  to  an  inextensible  microlayer  with  wave  reflection 
at  the  boundary  neglected.  However,  the  wind  input  inside  slicks  is 
required  before  a  realistic  representation  of  the  waves  inside  the  slick 
bands  can  be  modeled.  Investigation  of  this  wind  input  is  beyond  the  scope 
of  the  present  work.  A  realistic  representation  of  the  damping  must  also 
take  into  account  the  concentration.  Intuitively,  a  slick  should  approach 
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Figure  7. 
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The  Exponential  Decay  Time  for  Both  the  Free  Surface  and  the 
Inextensible  Microlayer  as  a  Function  of  Bragg  Wavelength 
Between  4  cm  and  24  cm. 
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an  inextensible  microlayer  only  near  it's  center  with  a  transition  to  a 
relatively  surfactant  free  surface  with  distance  away  from  the  center.  The 
microlayer  is  essentially  elastic  in  this  transition  region. 
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6.0  SUMMARY 


The  limited  measurements  reported  here  show  that  the  peaks  in 
fluorescence,  related  to  the  peak  concentration  of  Chlorophyll  A,  are  co¬ 
located  with  the  damping  of  the  short  gravity  waves  at  both  4  cm  and  24  cm 
wavelengths.  The  maximum  damping  ratio  between  the  24  cm  and  4  cm  waves  is 
at  least  16  for  a  wind  direction  of  about  30°  with  respect  to  the  radar 
look  direction  and  wind  speeds,  3-4  m/s. 

Imaging-radar  slick  signatures  displaying  narrow  bright  bands  on 
the  leading  edge  of  a  series  of  wider  dark  slick  bands  have  been 
interpreted  as  due  to  short  wave  straining  by  internal  waves.  We  assert 
that  the  wide  dark  bands  following  these  bright  bands,  when  observed,  are 
due  primarily  due  to  wave  damping.  That  is,  these  dark  bands  are  related 
to  concentration  of  the  marine  surface  microlayer  material.  Therefore,  the 
radar  roughness  of  slicks  are  related  to  both  short  wave  damping  and  wave 
straining  of  the  short  gravity  waves. 
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ABSTRACT 


We  study  the  feasibility  of  using  a  dual  linear-polarized  radiometer 
operating  at  6  GHz  in  conjunction  with  a  scatterometer  as  a  method  of  extracting 
sea  surface  temperature.  A  model  originally  developed  by  Wilheit  [1979]  is  used 
to  investigate  the  sensitivity  of  radiometric  temperature  to  changes  in  wind 
speed,  U,  at  a  height  of  20  m;  that  is,  the  accuracy  of  sea  temperature  depends 
directly  on  the  accuracy  of  radiometric  temperature.  These  changes  are  corrected 
using  this  extracted  wind  speed  based  on  a  scatterometer  wind  accuracy  of  ±  1 
m/s.  We  show  that  a  vertically  polarized  radiometer  is  less  sensitive  to  wind 
speed  changes  than  a  horizontally  polarized  one  near  an  incidence  angle  6  =  60°. 
Furthermore,  we  quantify  the  wind  variation;  ATr'^/AU  <  0.12°C/m/s  when  U  >  10  m/s 
and  ATb”/AU  <  2.0°C/m/a  when  U  >  6  m/s  for  vertical  and  horizontal  polarization, 
respectively.  Thus,  the  resulting  error  is  ±  0.12°C  and  ±  2.0°C  for  and  Tg”, 
respectively.  These  results  suggest  that  a  scatterometer  can  be  used  to  correct 
the  wind  speed  dependence  of  a  vertically  polarized  radiometer  operating  near  6 
GHz  using  a  pencil  antenna  beam  pointing  at  an  incidence  angle  near  57°. 
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1.0  INTRODUCTION 


Extraction  of  sea  surface  temperature  using  passive  microwave  radiometric 
measurements  in  conjunction  with  microwave  scatterometer  measurements  is  a 
technique  which  has  not  been  exploited.  The  technique  uses  the  scatterometer 
extracted  surface  wind  speed  to  correct  the  wind  dependence  of  the  radiometric 
brightness  temperature.  The  accuracy  achieved  depends  on  the  accuracy  of  the 
radiometric  brightness  temperature  measurement  itself,  the  accuracy  of  associated 
models,  and  the  accuracy  of  the  extracted  wind  speed  based  on  the  scatterometer 
measurements . 

The  radiometric  temperature  of  an  ocean  surface  depends  primarily  upon  its 
sea  surface  temperature,  salinity,  surface  wind  speed  (surface  roughness),  and 
foam  cover.  The  radiometric  temperature  of  a  smooth  surface  does  not  depend  on 
the  wind  speed  or  foaun  cover.  More  precisely,  foam  cover  is  not  important  for 
wind  speeds  less  than  about  7  m/s.  The  wind  speed  U  is  referenced  to  a  height  of 
20  m  since  most  surface  measurements  have  been  acquired  at  this  height.  The 
radiometric  temperature  of  a  smooth  sea  surface  depends  on  sea  surface 
temperature  and  salinity.  This  dependence  is  known  accurately  through  its 
relationship  with  the  complex  dielectric  constant  of  sea  water  [Klein  and  Swift, 
1977].  The  surface  wind  influences  the  radiometric  temperature  through  the  ocean 
wave  spectrum.  However,  measurements  are  not  fully  consistent  with  each  other  or 
with  the  theory  especially  when  foam  is  present  [Bollinger,  1971;  Monahan  and 
O'Muircheartaigh,  1980;  Smith,  1988]. 

This  study  is  limited  to  determining  an  optimum  radiometric  frequency  for 
active/passive  sensor  system  and  the  investigation  of  the  emissivity  dependence 
on  wind  speed.  The  model  examined  depends  on  atmospheric  md  ocean  parameters 
discussed  above.  However,  it  does  not  include  a  discussion  of  radiometric 
measurements  usually  taken  near  18  and  37  Ghz  to  correct  for  water  vapor  and 
clouds  which  are  most  important  at  aircraft  and  spacecraft  altitudes  [Guissard 
and  Sobieski,  1987]. 
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2.0  MODELS 


The  radiometric  temperature  Tg  for  a  near  surface  down-looking  radiometer 
over  a  rough  half  space  (wind  sea)  based  on  thermodynamic  eq[uilibrium  is  related 
to  the  sea  surface  temperature  Ts  (in  *K)  by. 


T,  =  e(S,aT,)T, 

where  S  is  the  salinity  of  the  ocean,  U  is  the  wind  speed  at  20  m,  e  is  the 
emissivity  of  the  rough  ocean,  and  Ts(°C)  =  Ts(°K)  -  273.16  .  The  emissivity  for 
a  smooth  sea  at  normal  incidence  is  given  by. 


=  1 


I 


/ 


(2) 


where  c  is  the  complex  dielectric  constant  of  sea  water. 

The  dielectric  constant  e  of  sea  water  is  dependent  on  the  ionic 
conductivity  a  in  mhos/m,  relaxation  constant  r  in  sec.,  and  static  dielectric 
constant  €,  (Debye,  1929]  as  given  by. 


e  =  ^ (J) 

1  +  {j  03  ry  '  ‘  (oCg 

where  £„  =  8.854  x  lo'^  is  the  permittivity  of  free  space  in  farads/m,  c.  is  the 
permittivity  at  a  high  frequency,  w  =  2irf  is  the  microwave  radian  frequency  with 
f  in  Hz,  a  is  an  empirical  constant  and  a,  r,  and  e  are  all  functions  of  S  and 
Tj.  Empirical  values  of  c,  r,  and  €  as  a  function  of  S  and  T*  were  given  by  Klein 
and  Swift  [1977]  for  the  microwave  frec[uencies  1.43  and  2.65  GHz.  We  use  these 
empirical  values  to  investigate  a  suitable  radiometric  frequency  for  which  the 
emissivity  dependence  on  salinity  is  negligible,  see  section  III.  The  selection 
of  this  frequency  allows  the  simplification  of  equation  (1)  resulting  in  T* 
depending  on  only  U  and  Tj, 


=  e(U.T,)T,. 

We  now,  in  principle,  can  extract  the  wind  speed  U  and  surface  temperate  Ts  using 
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equation  (4)  and 


where  is  the  scatterometer  backscattering  cross-aection  and  f '  is  the  inverse 
of  the  function  f  given  by, 

=  f(U)  =  aU» 

where  a  and  b  are  empirical  constants  dependent  on  the  scatterometer  frequency 
and  the  wind  direction.  The  extracted  winds  at  20  m  are  accurate  to  about  ±  1  m/s 
baaed  on  aircraft  circle  flights,  upwind,  crosswind  and  downwind;  see,  e.  g., 
Donelan  and  Pierson  [1987]. 

The  emisaivity  dependence  on  U  and  Ts,  given  in  equation  (4),  has  been 
modeled  quantitatively  by  Wilheit  [1979].  This  numerical  model  originally 
developed  by  Stogryn  [1967]  has  been  simplified  by  Wilheit  based  on  an  isotropic 
normal  distribution  for  the  ocean  wave  slopes.  Wilheit  obtained  good  agreement 
between  the  model  and  radiometer  measurements  based  on  this  simplification  with 
the  variance  of  the  slope  increasing  linearly  with  the  radiometric  frequency. 
Furthermore,  for  high  wind  speeds,  he  also  included  the  effect  of  foam  on  the 
radiometric  temperature.  He  treated  foam  as  partially  obscuring  the  surface 
independent  of  polarization.  Equation  (4)  and  this  model  are  used  to  estimate  the 
emissivity  dependence  on  U  and  Tj. 


87 


3.0  RESULTS 


The  optimum  radiometric  frequency  is  selected  such  that  both  the  emissivity 
dependence  on  salinity  is  negligible  and  the  atmospheric  effects  due  to  water 
vapor  and  clouds  are  minimized.  This  results  in  the  lowest  frequency  consistent 
with  negligible  dependence  on  salinity.  Blume  et  al.[1978]  used  radiometers 
operating  at  1.43  and  2.65  GHz  to  measure  salinity  with  aircraft  flights  over  an 
estuary.  Since  both  of  these  frequencies  are  sensitive  to  salinity,  an  optimum 
frequency  must  be  larger  than  2.65  GHz.  We  use  equations  (2)  and  (3)  and  Klein 
and  Swift's  results  in  the  range  3  <  f  <  10  GHz  to  investigate  the  minimum 
frecjuency  which  results  in  an  error  in  radiometric  temperature  of  less  than 
0.04°C  per  ppt  in  salinity.  The  salinity  of  surface  water  remains  relatively 
constant  for  long  periods  of  time  (near  34.6  ppt)  with  a  geographic  latitude 
variation  of  about  1  ppt  provided  specific  areas  are  excluded.  These  areas 
include  coastal  regions,  especially  estuaries  and  melting  ice  in  polar  regions 
[Williams,  1962]. 

The  results  of  our  investigation  are  shown  in  Figure  1  (a)  with  the 
emissivity  as  a  function  of  Tj  with  an  optimum  frequency  of  6  GHz;  while.  Figure 
1  (b)  illustrates  the  same  variation  at  2.65  GHz  for  comparison.  The  range  of 
salinity  was  varied  from  20  ppt  to  40  ppt  in  intervals  of  2  ppt.  Each  ppt 
interval  change  in  Figure  1  (a)  corresponds  to  <  0.00013  change  in  emissivity 
or  equivalently  <  0. 038'’C/ppt .  Each  ppt  interval  in  Figure  1  (b)  corresponds  to 
<  0.0007  or  <  0.22‘’C/ppt.  These  results  are  based  on  Tj  S  20“C.  Therefore  the 
optimum  radiometric  frequency  is  near  6  GHz  as  given  in  Table  2(a).  Indeed,  a 
typical  salinity  change  of  1  ppt  with  latitude  causes  a  change  in  Tj  of  O.OSO'C 
at  C-band  (6  GHz)  and  0.22'’C  at  S-band  (2.65  GHz)  at  normal  incidence. 

We  now  consider  the  emissivity  dependence  on  only  U  and  Tj  using  equation 
(4)  and  Wilheit ' s  model.  The  variation  of  the  emissivity  with  incidence  angle  for 
three  wind  speeds,  U  =  6,  10  and  20  m/s  are  given  in  Figure  2  (a).  Wind  speeds 
less  than  6  m/s  have  negligible  foam  cover.  The  upper  three  curves  correspond  to 
vertical  polarization  and  the  lower  three  curves  to  horizontal  polarization.  The 
vertically  polarized  curves  tend  to  intersect  at  an  incidence  angle  6  between  50° 
and  60°.  Therefore  these  incidence  angles  correspond  to  a  minimum  dependence  on 
wind  speed  which  is  the  same  result  as  previously  reported  by  Stogryn  [1967].  The 
percent  of  effective  foam  cover  for  wind  speeds  from  8  to  20  m/s  based  on  the 
model  are  given  in  Table  1.  Figure  2  (b)  shows  the  same  results  as  Figure  2  (a) 
except  the  differential  emissivity  is  displayed  on  the  vertical  axis.  This 
normalized  differential  emissivity  6e/e,  in  percent  is  defined  as. 
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TABLE  2(a). 

Radiometric  Brightness  Temperature  Sensitivity  to 

Salinity  at  Normal  Incidence 

OPERATING  FREQUENCY  I 

C-BAND  (6  GHz) 

S-BAND  (2.65  GHz) 

ATb/AS 

0.038‘’C/ppt 

0.22‘’C/ppt 

TABLE  2(b). 

Radiometric  Brightness  Temperature  sensitivity  to  Wind  Speed, 

S  =  34  ppt  and  6  =  60°. 

C-BAND  ( 6  GHz ) 

POLARIZATION  I 

VERTICAL 

HORIZONTAL 

ATb/AU 

<  0.12'’C/m/8 

<  2.0“C/m/8 
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TABLE  1 


WIND  SPEED  AND  ASSOCIATED  PERCENT  EFFECTIVE  FOAM  COVER 
BASED  ON  WILHEIT  [1979] 


U(m/B) 

FOAM ( % ) 

6 

0 

8 

.33 

10 

.99 

12 

1.65 

14 

2.31 

16 

2.97 

18 

3.63 

20 

4.31 
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6e(e) 

eJB) 


(es(d)  -  e(d)) 

es(e) 


% 


(7) 


where  e,(6)  is  the  emissivity  of  the  smooth  sea  as  a  function  of  incidence  angle 
and  e,„{0)  and  e,v(0)  are  the  corresponding  vertical  and  horizontal  polarizations 
of  e,(0)  given  by. 


/ 


cos6  -  -  sir^6 

cos  6  +  ye  -  sir^d 


/ 


2 


(8) 


=  J 


/ 


e  cosd  -  /e  -  sin^d 
e  cosd  +  ye  -  sin^d 


I 


(9) 


The  three  lower  curves  correspond  to  vertical  polarization  and  the  three  upper 
curves  correspond  to  horizontal  polarization.  The  differential  emissivity  for  the 
vertical  polarization  consistently  shows  a  lower  percent  variation  as  a  function 
of  incidence  angle  and  wind  speed. 

He  next  consider  the  emissivity  dependence  on  sea  surface  temperature  for 
a  range  of  wind  speeds,  U  6  -  20  m/s  at  0  =  60°.  Figure  3(a)  shows  the 
emissivity  dependence  on  U  and  T^  for  horizontal  polarization  whereas  Figure  3 
(b)  shows  this  same  dependence  for  vertical  polarization.  The  primary  result  is 
that  the  vertical  polarized  Tg''  is  only  weakly  dependent  on  wind  speed  when  6 
m/s  <  U  <  10  m/s  whereas  at  higher  wind  speeds  the  change  in  emissivity  is  Ae/AU 
<  0.00042 /m/s  or  in  terms  of  7^  is  ATg'^/AU  <  0.12°C/in/s.  In  contrast,  the  change 
in  emissivity  for  horizontal  polarization  is  Ae/AU  <  0.007/m/s  or  ATg  /AU  < 
2.0°C/m/s  when  U  >  6  m/s.  These  Tg  sensitivities  to  wind  speed  are  summarized  in 
Table  2(b).  The  maximum  error  is  ±  0.12°C  and  ±  2.0°C  for  Tg^  and  Tg" 
respectively  based  on  the  scatterometer  wind  accuracy  of  ±  1  m/s  assuming  no 
errors  in  the  model  and  the  measurement  of  the  Tg^  and  Tg". 
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Emissivity , 


horizontal  polarization 


Figure  3.  EmisBivity  Variation  with  Sea  Surface  Temperature  T*  Parameterized 
in  Surface  wind  Speed,  U  =  6  -  20  m/s,  6  «  60*  and  Salinity,  S  =  34 
ppt;  (a)  Horizontal  Polarization,  and  (b)  Vertical  Polarization. 


4.0  SUMMARY 


The  model  results  suggest  an  experiment  using  a  tower-based  radiometer 
operating  at  6  GHz  with  measurements  taken  simultaneous  at  horizontal  and 
vertical  polarizations  using  a  pencil  beam  antenna  selected  at  an  incidence  angle 
near  57°.  These  measurements  would  be  accjuired  in  conjunction  with  scatterometer 
rind  wind  anemometer  measurements  in  an  area  where  high  winds  and  ocean 
temperature  fronts  frequently  occur  allowing  the  validation  of  the  models  and  the 
measuring  technique. 
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